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Abstract. 


Snowmelt modelling in different climatic regions and various types 
of terrain is of major importance in snow hydrology. 


During 1979-82 a field research programme was carried out at 
Dyrdalen (60°21'N, 5°31'E), an area of rough terrain in the mari- 
time climatic region near the coast of western Norway. 


Runoff from 9 mé 


four locations, two of which were situated above the timberline, 
at 437 and 632 m a.s.1., respectively. The other two stations 
(both at 285 m a.s.1.) were located in an area partly timbered by 
deciduous trees. Wind speed, air temperature, air humidity, 


plots of snow cover was recorded by lysimetry at 


radiation, and precipitation data were also recorded. 


When the melting rate (= snowpack runoff-rainfall) and the net 
radiation were measured, the turbulent heat exchange between the 
snowpack and the atmosphere was computed as a residual from the 
energy balance equation of the snowpack. These computed values 
were used to find "optimal" empirical constants in aerodynamical 
equations expressing the turbulent fluxes as functions of the 
wind speed and the temperature/vapour pressure differences between 
the measurements 2 m above the ground (0.6 - 1.8 m above the snow 
surface) and the values at the surface. These empirical constants 
agree reasonably well with those found by other investigators. 
The obtained constants proved to be near optimal for the area 
above as well as below the timberline. 


This energy balance model yielded higher model efficiency than a 
degree-day model for locations above the timberline, while the 
degree-day model was the better in the wooded area. The degree-day 
model, however, has to be tuned to the actual type of vegetation. 


The efficiency of the models proved to be high also when data for 
wind speed, air temperature, and air humidity were extrapolated 
from records at stations up to a distance of some 10-20 km. High 
quality estimates of net radiation are needed to get reliable 
results during fine weather, while rather rough estimates are 
sufficient during cloudy weather. 


List of the most frequently used symbols. 


a = albedo (%, fraction of unity) 
C = cloud cover (%, fraction of unity, octants) 
Cy = bulk transfer coefficient 
Cy = spesific heat of dry air (J kg ko!) 
.* Pa a 
" _ (mb K — ) 
0.622 L. 
e = vapour pressure (mb ) 
€, = vapour saturation pressure at the surface (mb ) 
€ = emissitivity of snow (fraction of unity) 
g = acceleration of gravity (ms~2) 
k = degree-day factor (mm day 'ec7!) 
ky = ratio between ground level records of 
precipitation and records some height 
(1.6, 2 m) above the surface) 
k = von Karman's constant (0.41) 
L, = latent heat of sublimation (J kg!) 
n = number of observations 
P, = air pressure (mb ) 
P = true precipitation (mm) 
P, = observed precipitation by (mm ) 
Belfort instrument 
Py = observed rain amount by lysimeter (mm) 
P., = observed precipitation (mm) 
Po = observed rain amount by OTA (mm) 


ground level instrument 


Q, = latent heat flux (wm?) 
Qa = extraterrestrial global radiation (wm 2) 
Q, = sensible heat flux (wm?) 
Q,, = incoming longwave radiation (um) 
Qn = computed energy consumption for snowmelt (wm?) 
Qy' = observed a it a (wn?) 
Qy = net radiation (Wwm-*) 
Qo = global radiation (wm 2) 
ig = bulk Richardson number (dimensionless) 
Re = critical Richardson number 
R2 = model efficiency . 
p, = density of air (kg m7 3) 
p,, = density of water (kg m9) 
SD = standard deviation 
Si = predicted daily snowmelt by the degree- (mm day!) 
day model 
o = Stephan Bolzmann constant (wm 2x74) 
Oo, = residual error 
T = air temperature ("Gs #2 
t = time of exposure to open air of the snow (days) 
at the surface 
T, = threshold temperature LG) 
us = ~wind speed (ms~!) 
Z) = surface roughness (cm) 
Z = recording level (1) 


Subscripts a: air level 0 : surface 


1. Introduction. 


Snow and ice, which are common features of the Norwegian environ- 
ment, play a significant role in the hydrological cycle. The 
accumulation and melting of snow strongly affect runoff from high 
altitude and/or high latitude catchments. An adequate understand- 
ing of the physics of snowmelt is therefore of vital importance 

to water resources management in such areas. High precipitation 
amounts, with significant snow accumulation and frequent snowmelt 
caused by extremely diverse weather conditions, make the western 
part of southern Norway well-suited for studies of certain relation- 
ships between snowmelt and weather. Knowledge of the snow distribu- 
tion, snowmelt and water movements in snow are essential in plan- 
ning hydroelectric power production and for predicting snowmelt 
floods. 


During winter, snow is the most frequent type of precipitation in 
extensive parts of the country, and in the mountains the snow 
accumulation season usually lasts from October to May. The snow- 
melt usually takes place in March and April in the lowlands, and 
in May and June in the mountains. The times are somewhat later 

at high latitudes. Near the western coast the snow cover is highly 
unstable due to the maritime location. Snow accumulation events 
and snowmelt periods are both common throughout the winter. At 

low altitudes (<300 m a.s.1.) near the sea, winter precipitation 

is generally rain or sleet; at medium altitudes (300-600 m a.s.1.) 
there is usually a main snowmelt period in spring, and several 
periods with heavy rain or snowmelt during the snow accumulation 


season. 


The snow is distributed in variable patterns depending mainly on 
terrain parameters and surface wind. Large variations in snow 
depth are found on horizontal scales of 10-1000 m. Newly fallen 
snow has low density, i.e. 50-300 kg-m 3, the highest values 
being found for wet or wind-packed snow. The density increases 
with time to 300-600 kg m3, ag the snow crystals grow. The time 


scale of the density increase depends mainly on the level and 


variations of the air temperature. Heavy rain episodes may 
accelerate the process. A wet snow pack (0°C) holds a certain 

amount of liquid water against the gravitational force. This 
amount, expressed as a percentage of the total snow weight, is 
called the free water of the snow pack. The maximum free water 
content is called the liquid water capacity. U.S.A.C.E. (1956) 
report it to be 2-5% by weight. Colbeck (1974) and Nyberg and 
Harsmar (1971) give the value 4%. The actual value obviously 

varies considerably. In addition, during rain and snowmelt episodes 
the snow pack contains liquid water percolating downward. 


A dense and grainy snow pack with free water at liquid capacity 
is said to be ripe. Any rain or melt water entering the snow 
surface now leaves the bottom as discharge. The percolating time 
varies with depth and nature of the snow pack, but it is usually 
on the order of a few hours or less (Colbeck, 1972). 


The melting process requires energy which is supplied as net 
radiation, and sensible and latent heat from the atmosphere. 

The ground heat flux and heat from rain also make minor contribu- 
tions. The air temperature may be a useful index of snowmelt 
Since it is often correlated with net radiation, and with the 
turbulent fluxes of sensible and latent heat. 


Accordingly, meteorology plays a major role in snow hydrology. 
The subject is also closely related to glaciology. U.S.A.C.E. 

in their work Snow hydrology (1956) discuss the subject in its 
entirety. Extensive studies of water movements in snow have been 
made in the U.S.A. and Canada during the last 20 years (Colbeck, 
1972, Wankiewicz, 1978). The problem of snow distribution in 
Norway has been discussed by Tveit (1979), Andersen and Odegaard, 
(1980) and others. Operative models for water field discharge, 
including snowmelt, have been developed in several countries. In 
Scandinavia, Bergstrom (1976) developed a model for snowmelt 
using a temperature index method. This model is the most frequently 


used in Scandinavia today. 


Sverdrup (1936) was the first to carry out extensive studies of 
the energy exchange between a snow surface and the atmosphere. 
Wallén (1948) examined the energy balance of a Swedish glacier 
surface during five summer seasons and found the relative contri- 
butions of net radiation, sensible heat and latent heat to be 
55%, 29%, and 16%, respectively. Paterson (1969) presented a 
table covering many glacier examinations which showed highly 
variable contributions. U.S.A.C.E. (1955, 1956) carried out 
lysimeter studies of snowmelt in order to develop a model based 
on the energy balance of the snow cover. Liljequist (1956-57) 
made extensive studies of radiation and turbulent fluxes at an 
Antarctic field. Studies of this problem were also carried out 

in the U.S.S.R., published by Kuzmin (1961, in English 1972). 
Krauss (1966) examined glacier ablation and produced curves of 
ablation for different values of air pressure, wind speed, air 
temperature, air humidity, global radiation, albedo, and surface 
roughness. Holmgren (1971) studied the energy balance of an Arctic 
glacier during three summer seasons; he found the most favourable 
conditions for snowmelt here to be strong wind and advection of 
warm and humid air. Anderson (1964, 1967, 1968, 1973, 1976) 
carried out further studies to solve the problem of complete 
energy balance. He also took into account stability corrections 
for the turbulent heat fluxes, and modelled the energy transfer 
within the snowpack.. In Norway snowpack discharge has been 
recorded by lysimeter (Tveit, 1977), and changes of snow weight 
by snow pillow (Furmyr, 1975), but the data have not been used 
for energy balance studies. Snowpack discharge may also be calcu- 
lated from regular observations of snow depth and density. For 
areal estimates, field discharge may also be directly recorded. 


The question of whether an energy balance model can be an improve- 
ment on existing snowmelt models is of considerable interest. 
Therefore, a field research programme supported by the Norwegian 
National Committee for Hydrology and University of Bergen was 
carried out near the coast of western Norway. In this paper the 


energy balance model and a temperature index model are tested 
against lysimeter data at a woodless index point during periods 
where the snow pack could be considered ripe and isothermal, 
except for daily variations. It is also attempted to assay how 
snowmelt is modified in areas partly timbered by birch and other 
deciduous trees. The data requirement and tuning of the models 
are also discussed. Preliminary results on some of these topics 


are given by Harstveit (1981). 


2. Site Description and Climate. 
2.1. Site Description 


During 1978-1979 two stations, Dyrdalsvatn (437 m a.s.1.) and 
Austlihylla (632 ma.s.1.), were established in the Dyrdalen 
area (Fig. 2.2) forrecording snowmelt. 


Dyrdalen (60°21'N, 5°31'E) is situated above the timberline, 
about 11 km ESE of Bergen. The precipitation catchment area 
(3.34 km*, 435-806 ma.s.1.) is a rough mountainous area (Fig. 2.3), 
but on a smaller scale (%100 m) the terrain is smoother. The 
main direction of the valley is SSW-NNE, and it is situated 


above the timberline. 


The soil and morain masses extend to a depth of 0-2 m. Boggy 

ground and low vegetation of grass, heather, and moss are pre- 
dominant. A more detailed description of the field in general, 

the geology, and the vegetation is given by Gjessing et al. (1980). 


During the winter 1980-81 the station network was extended with 
two stations at Frotveit (285 m a.s.1., Fig. 2.2). These stations 
are situated about 5 km SW of Dyrdalen, and the distance between 
them is 70 m. Though the terrain in this area is less rough, 
there are considerable variations on a smaller scale (Fig. 2.5). 
The direction of the valley is SW-NE, and it is partly timbered 
by birch and other deciduous trees, and some small area with 
densely planted spruce. The terrain and vegetation around the 
stations are illustrated in Fig. 2.5. One station is located at 
an open place, and the other one in a grove consisting mainly of 


birch and some older and rowan. 


Data from weather stations (1-3, Fig. 2.1) operated by the Norwegian 
Meteorological Institute are also used in this investigation, 
together with radiation data from Bergen - Florida, published by 
Geophysical Institute, University of Bergen (Radiation Yearbook 
1-18). 


Figs Zals Map of the area. 


1 Bergen, 2 Flesland, 3 Kvamskogen, 4 Austlihylia, 


5 Dyrdalsvatn, 6 Frotveit. 


DYRDALEN 


2.3. View of Dyrdalen. The snow amounts in March 1979 (below). 
1. Dyrdalsvatn (Weather station) 
2. Dyrdalsvatn (Lysimeter) 
3. Austlihylla 


Fig. 2.4. Austlihylla snowmelt station on June 13, 1979. 
The snow-drift is illustrated. (The Belfort precipitation 
gauge which had been damaged by snow weight during this 


winter was moved to the right of the lysimeter in July.) 
1. Weather station 


2. Lysimeter 
3. Belfort rain gauge 


2s Bs 


The woodless (above) and wooded (below) site at 


Frotveit. 


1, Stand Tor 2. Snow lysimeter 3. Rain gauge 
Radiation records 
Anemometer 
Hygrometer 
Thermometer 


2.2. Climatic Review. 


Western Norway is exposed to prevailing westerly and south- 
westerly winds which bring mild and moist air from the North 
Atlantic Ocean over the country all year. The standard normal 

of the annual air temperature is about 7 to 8°C in the lower 
districts, while the monthly values are about 0°C in January 

and 15°C in July. In March the standard normal is about 2 to 3°C, 
5 to 6°C in April, and 10°C in May. The standard deviation of 

the monthly mean temperature is 2 to 3°C in winter, and 1 to 1.5°C 
during the rest of the year; it is highest for locations exposed 
to cold drainage, and lowest near the sea. The annual and monthly 
values for Bergen and Kvamskogen are illustrated in Fig. 2.6. The 
temperatures at Dyrdalen should be close to those of Kvamskogen. 


The precipitation is high throughout the year (Fig. 2.7). 34% of 
the annual amount falls in September-November, and 17% in March- 
May. There is a maximum zone of precipitation in the mountains 
30-60 km inland (Forland, 1980 a). The local variations, however, 
are large due to the rough terrain. The standard normal of the 
annual precipitation, which is 1958 mm in Bergen and 2755 mm at 
Kvamskogen, is estimated to 3300 mm at Dyrdalen (Abildsnes, 1980). 
The highest precipitation amounts recorded at Kvamskogen are: 

151 mm in 24-hours (07 p.m. - 07 p.m.), 864 mm over one calendar 
month, and 4506 mm over one calendar year. 


It is difficult to give an adequate description of the wind 
conditions in the area due to the rough terrain. Fig. 2.8 shows 
wind roses from Flesland, which is a relatively exposed wind 
station, near the coast. The monsoon effects are significant 

with predominating northwesterly to northerly winds in summer 

and easterly to southeasterly winds in winter. The frequent 
southeasterly to southerly winds throughout the whole year reflect 
the cyclones coming from southwest. The surface wind speed at 
Dyrdalsvatn should be of the same strength in such situations, 


For winds from east to southeast, Dyrdalsvatn is sheltered by 
high mountains. Austlihylla is exposed to these directions, and 
also to northwesterly winds. Frotveit is far more sheltered than 
the two higher stations, being most exposed to southwest and 


northeast. 


The cloud cover in Bergen was in average (1901-60) 5.6 octants. 
The annual distribution (Fig. 2.9) shows a minimum in May. The 
relative sunshine duration shows a maximum i April-May and a 
minimum in December. The average duration during 1965-79 was 
31% of the possible value. 
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Mean air temperature 1931-60 
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Monthly and annual means of precipitation at Bergen- 
Fredriksberg. and Kvamskogen during the period 1931-60. 
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Fig. 2.8. Mean seasonal and annual wind roses for Flesland, 1956-65. 
% of observations in each of 12 direction sectors; rela- 
tive fraction of observations in each of 3 speed 
(Beaufort, 8) intervals; % calm at center. Summer: 
June-August. Winter: December-February. 


+——+ Relative duration of sunshine,-= 
= Ovevece o Cloud cover,C (Olseth, 1980) C 
0.4 
t—t—) ae 
J NH 
son 
0.3 \. 
~ 
nN 6.0 
2 *\.0 
0. B79, 991710 0,0, 00 Ors 
° 2 + 
nia 5.5 
0.1 ‘ . fe) | 
ad 5.0 
J F M AM JS J A S O N ODO 


Fig. 2.9. Monthly values of the cloud cover (1901-60) in octants, 
and the relative duration of sunshine (1965-79). 
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3. Instrumentation and Methods. 


3.1. Instruments, Data and Data Quality. 


In the Dyrdalen-Frotveit area records were made at four stations 
(Table 3.1). All sensors were placed 2 m above snow-free ground, 
except for the precipitation gauges (1.6 m), and the snow lysi- 

meters (ground level). The micrometeorological instruments had a 
mean height of roughly 1.3 m above the snow surface during snow- 


melt. 


The data were recorded on magnetic tape every 30 minute by data- 
loggers manufactured by Aanderaa Instruments, Bergen. 


Table 3.1. 
Instrumentation and recording periods during snowmelt. 


Instrument Austlihylla Frotveit I Frotveit II Dyrdalsvatn 
Snow lysimeter 1980-82 1981-82 1981-82 1979-82 
Precipitation gauge 1980-82 1981-82 1979=82 
Thermometer 1979-82 1981-82 1981-82 1979-82 
Hair hygrometer 1980-82 1981-82 1981 1979-32 
Anemometer 1980-82 1981-82 1981-82 1979-82 
Net radiometer (Total radiation) 1981 198] 

Net radiometer* 1980-82 
Pyranometer (Upfacing) 1981 1981 1979-82 
Pyranometer (Downfacing) 198] 1979-82 


* Adaptor covering lower hemisphere 
I: woodless site 
II : wooded site 


The snow lysimeter 


Snowpack discharge was recorded by a snow lysimeter, described 
by Tveit (1977) and illustrated in Fig. 3.1. Two polyester vats 
with walls 0.15 m high were placed side by side at ground level 


covering a total area of 9 m2, The vats were allowed to snow down, 


INSULATED COVER WATER LEVEL RECORDER 


COLLECTION VATS | 
FOR PRECIPITATION | 


DRAINPIPE 
DISCHARGE COLLECTOR | 


{ BAILING SEESAW 
DRAIN 


Fig. 3.1. Principle of the snow lysimeter used at Dyrdalen 
and Frotveit 1979-1982. 


and during snowmelt the discharge was recorded by measuring the 
water level of a tank of capacity 112.5 1 (12.5mm discharge). This 
collector tank is automatically emptied by a bailing seasaw when 
the water level reaches the top. It is placed well beneath ground 
level and insulated against frost. To avoid boundary effects, 
the data series were truncated either when the snow coverage in 
the lysimeter fell to less than 90% of the bottom area, or when 
the corresponding water equivalent fell to less than 50 mm. 


The snow cover near to the Dyrdalsvatn snow lysimeter was almost 


complete during all the snowmelt periods recorded. 


- 22 - 


Due to percolation resistance in the snow there is a time lag 
between actual and recorded snowmelt. For daily totals during 
rain-free snowmelt episodes, however, the well-defined diurnal 
wave makes time lag corrections unessential, except at Austlihylla 
where a deep snow cover prevailed. During rain-on-snow events 

some correction of the timing of the discahrge was necessary, 
mainly at Austlihylla, and was carried out by comparing lysimeter 


discharge and rain gauge records. 


U.S.A.C.E. (1955) report that there may be a net flow of liquid 
water over the walls of a lysimeter. When that flow was eliminated 
(using a technique of adapting the lysimeter edges to the snow 


’ 


depth), their data were very accurate. 


Some checks of lysimeter data at Dyrdalsvatn are summed up in 
Table 3.2. The water equivalent, w.e., determined by measuring 
snow depth and density on the lysimeter vats, is added to recorded 
precipitation, Ps » and estimated condensation, CD , accumulated 
until the lysimeter is completely empty. When comparing this sum 
to the accumulated lysimeter discharge for the same period, only 
small deviations are found, indicating that flow across the lysi- 


meter edges at Dyrdalsvatn is a rather small source of error. 


By a rough guess the standard error of daily snowmelt is estimated 

to be within 5%. The mean free water content of the Snowpack is 
Supposed to be 4% (Capt. 1). In some cases the snow cover is strongly 
layered, and the free water content may deviate considerably from 
this mean value (Colbeck, 1972). The error will then be larger. 


At Austlihylla the terrain is sloping (Fig. 2.2), and the snow 
cover is deep due to the local snow-drifting, conditions which 
make a significant flow across the lysimeter more likely. However, 
due to breaks in the lysimeter data series, there is no avaiable 
snowmelt period for checking at this station. 
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Table 3.2. Results from checks of lysimeter data of Dyrdalsvatn (see text). 
Numbers in brackets are rough estimates of the standard errors. 


Date of measuring, w.e. 
Parameter May 10, 1979 April 16, 1980 April 19, 1982 


Number of days till 


the lysimeter was dry \e a : 


Water equivalent 


fii ®a ts. tith 440 (18) 273 (14) 98 (4) 
Precipitation 210 (20) 56 (6) 28 (3) 
(P). mm 

Estimated condensation 10 (4) -6 (2) 0 (1) 
(CD), mm 

w.e. + P a + CD 660 (27) 323 {13} 126 (5) 
Discharge (Q), mm 650 (13) 295 (6) 129 (3) 


Table 3.3. Comparison of the daily discharge records from separated 
lysimeter vats (a,b) at Frotveit during snowmelt 1982. 
In brackets standard deviation of daily values. 


Woodless site Wooded site 
Number of days 9 11 
Daily discharge, Q, (mm) 14.8 (12.1) 16.3 (13.8) 
Daily discharge, Q, (mm) 1633 (1167) 18.4 (12.9) 


Correlation coefficient Fat. = 0.998 Pep = 0.999 
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In 1982 the lysimeter vats at both stations at Frotveit were 
equipped with separate recorders. Results from the available 

data are shown in Table 3.3. The correlation coefficients between 
records for adjacent vats are very high. The systematic deviations 
are small, especially for the open site. For the wooded site this 
deviation is probably caused by differential shading by surrounding 
trees. The analysis indicates that the quality of the discharge 
data from the lysimeters at Frotveit is high. The standard error 
should be as low as for Dyrdalsvatn, suggested to be within 5%. 


Rain records. 


Records of liquid precipitation generally include several sources 

of error. A comprehensive discussion is given by Dahlstrom (1970), 
and Killingtveit (1976). The aerodynamical effect usually repre- 
sents the largest error. Local effects (vegetation, wind conditions) 
at the station may make the records unrepresentative of the 
surrounding area. Other error sources are splash out or into the 
collector, wetting of the collector, evaporation loss and general 
interpretation errors introduced when deciphering the data. 


Precipitation records were made by weighing pluviographs, type 
Belfort equipped with Alter type wind shields, 1.6 m above snow- 
free ground at Dyrdalsvatn, Austlihylla and Frotveit (open site). 
At Dyrdalsvatn precipitation was also recorded by a tipping bucket 
pluviograph (OTA) at ground level during parts of the summer 
seasons. The lysimeters also functioned as ground level rain 
recorders during parts of the snow-free seasons. 


The ground level records are assumed to be free of aerodynamical 
effects. The integrated error of the OTA record should therefore 

be very small, and clearly within 0.5-1 mm for daily totals 
(Forland, 1980). The errors in the lysimeter rain data are within 

= 9 except for splasheffects (Killingtveit, 1976). For the Belfort 
records the error due to aerodynamical effects usually is consider- 
ably larger than the integrated error due to other effects es 1 mm 
for daily totals) (Fgorland, 1980). At Austlihylla only a few 


< OF 


records are available due to mechanical problems, so no analysis 
was carried out using those data. The few data indicate, however, 
that the aerodynamical effect is larger and more variable than 


at Dyrdalsvatn. 


Fig. 3.2 shows some results from an analysis of data recorded 
during the snow-free seasons of 1979 and 1980. Daily precipitation 
records, ee > 10 mm are grouped according to the amount, 3 in 

each group. The group mean values of Po Pp » andP, (defined inFig.3.2)ari 
calculated, and the ratios between corresponding group mean values 
are plotted against daily precipitation. Lines of linear regression 
are drawn, though neither of the regression coefficients are 
statistically significant. However, splash loss from the lysi- 
meter during heavy rain €an explain the negative slope of the 

lines 2 and 3. The figure indicates that the splash loss increases 
some 1-2% for each 10 mm increase in daily precipitation for totals 
between 10 and 60 mm. For small daily totals the splash loss is 
assumed to be zero due to the edges of the lysimeter which catch 
splash from light precipitation. Accordingly, for daily totals 
below 10 mm no splash loss from the lysimeter is assumed,. and 
equations 2 and 3 are not valid. For a daily total of 10 mm the 

OTA and Dyrdalsvatn lysimeter records only deviate with 0.4%. 

The splash loss is assumed to be zero during rain on snow events, 


and the equation 


P~P,.~% PP =KP (3.1.13 


LD 0) p B . 


where P is true precipitation, and K, is the ratio between 
ground level and 1.6 m level precipitation records, should then 


be valid. 


P 
Pron Fig. 3,2 mate = 1.09 and 1.06 (for Pe = 10 mm) for the two 

PB 
data sets, comparing data recorded at Dyrdalsvatn by the OTA gauge 
(P,) and the Belfort gauge (Pg). The difference may be due to 
different wind conditions during the days involved. Furthermore, 


Pi p/Pp = 1.10 for the set of highest = value. The average value 
B 


F4g: 


Lysimeter at Dyrdalsvatn 


1.086 +0.0005 P, ,r=0.19 


= 1.018-0.0014 B, ,r=-0.66 


Ratio | 
= 1.109-0.0013 Py ,r=-0.44 
1.20 | 
x 1 
—_-x-—OO 
1.10 x 
a 
3 
1.00 
a 2 | 
0.90 [ll eee 
10 20 30 40 50 60 = mm day? P, 
Lysimeter at Austlihylla 
= 1.064 ; Fr =-0:01 


1.039-0.0016 Py, r =-0.31 


= 1.110-0.0018 Py. r =-0.30 


3.2. Relation between ratio of recorded rain amount 
different recorders and daily totals of precipi 
Each data point represents a group mean value o 
observations. 


Total of precipitation recorded by OTA pluviograph, 
us " by Belfort pluviograph, 


by lysimeter, 


by lysimeter, 


60 70 80 mm day Py | 


between 
tation. 
f 3 


Dyrdalsvatn. 
Dyrdalsvatn. 
Austlihylla. 
Dyrdalsvatn. 
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of kK therefore is taken to be 1.08. Abildsnes (1980) found the 
value 1.08 between ground level and wind-shielded 2 m level 
precipitation records at Dyrdalsvatn. Furmyr (1975) found 1.09 
at Fillefjell. 

PLA ‘ 
For Austlihylla —— = 1.02 (for p, = 10 mm), suggesting a 2% increase 
in precipitation Senate from Dyrdalsvatn to Austlihylla. Therefore, 
Ky is put to 1.10, where Ky is analogous to Ky with the effects 
of the intersite distance included. When the data published by 
Abildsnes (1980, pp. 27) are used the precipitation amounts at 
Dyrdalsvatn and Austlihylla are found to be very equal on the average. 
There were, however, significant differences between individual 
episoages depending mainly on the wind direction. From Fig. 3.2. it 


appears that the dispersion of the a values are larger for the 
lysimeter at Austlihylla than for E Dyrdalsvatn, which should 


be expected since Pe is recorded at Dyrdalsvatn. The standard 
ceviations of oe for daily totals were respectively 17% (n=24) 
and 6% (n=30) of average ratios at Austlihylla and Dyrdalsvatn. 


During the snow-free and leafless season at Frotveit, daily precipi- 
tation above 10 mm was recorded for 10 days. The average daily 
total was only 16 mm and the standard deviation 4 mm for those 
days. Splash effects should be rather small due to tne lack o 

heavy precipitation records. Therefore the average value of _ , 
where Py is the lysimeter total and Pe is the total recorded 

by the wind-shielded Belfort rain gauge situated at tne woodless 
site, can be used directly as an estimate for the correction factor 
of Pe The values 1.04 and 0.98 are found for the open and wooded 
sites, respectively, with standard deviation of 5% and, 11%. The 
low factor at the open site should be due to the lower wind speeds 
in the area compared to Dyrdalen. For the wooded site the even 
lower factor 0.98 reflects the fact that some of the precipitation 
is intercepted on the trees. The standard deviation of 11% reflects 


that the interception is variable. 


To sum up, the daily totals recorded by the wind-shielded Belfort 
rain gauge at Dyrdalsvatn should be multiplied by 1.08 to estimate 
true precipitation at this site, and by 1.10 to give estimates 
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for Austlihylla. At Frotveit the totals recorded by the Belfort 
gauge should be multiplied by 1.04 to estimate the true precipita- 
tion at the woodless site, and by 0.98 to give estimates for the 
wooded site. These factors are averaged values. The actual values 
vary with the wind conditions, and for the wooded site, the inter- 
ception conditions.’Inpercents the standard deviations for the 
factor Ky (6%, 17%, 5%, and 11%) should be. reasonable esti- 
mates for the error due to the aerodynamical effect when a first 
order correction is made. The standard error, AP (mm) of "true" 


precipitation is then given by, 
é dL 
+ (PBR) 2 Umm), , (3.1.2) 


where APap (= 1 mm) is the integrated error due to other effects 
of the precipitation totals recorded at the Belfort gauge, Pp 


For heavy precipitation the second term under the radical is small 
AK 
(mm) or — 100 (%) 


p Kp 


and the standard error is given by, Kp . Pp 
K 


Precipitation during snowmelt. 


The precipitation form is related to the air temperature. In Fig. 3.3 
is shown the probability for liquid precipitation as a function of 
the air temperature, T , based on a study in the U.S.A. iY. Ss Ac. Eas 
1956). At T = 1.5°C the probability is 50%. Furmyr (1975) and 
Killingtveit (1976) both report 1.1°C as the temperature where the 
precipitation changes from snow to rain. 


Accordingly, the precipitation in the Dyrdalen-Frotveit area is assumed 
to be liquid when the air temperature exceeds 1.1°C. Additional 
information such as direct observation or obviously snow-covered 
pyranometers may in some cases improve this classification. Periods 
with large fluctuations in the precipitation form are excluded from 
the analysis due-to the errors introduced by estimating the rainfall 


contribution. 


R (%) | 

100 @ 

80 | 
Snow 

60 Rain 
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Fig. 3.3. Precipitation form in per cent as group mean values of 
precipitation within groupes of 0.56°C between -1.4 and 
-4.1°C (Data from Table 3 - 02.03, U.S.A.C.E., 1956). 


During precipitation the snowmelt rate is given by subtracting the 
rain amounts from the discharge. Standard errors of precipitation 
records then increase the standard errors of the snowmelt rate. 
The following example illustrates this effect.(for Dyrdalsvatn): 


The standard error of the discharge is assumed to be + 5%. For a 
daily snowmelt total of 40 mm this error makes the error of the 
snowmelt records to ~ 2 mm. If 40 mm rain falls in addition to the 
snowmelt assumed, the discharge error is 2 mm, and eq. (3.1.2) 


gives AP = 2.6 mm. The error of the snowmelt total then is 


+ ~ 
increased to ((2.6)°+42)2mm % 5 mm or 12%. 


For precipitation at air temperatures between O°C and 3°C, the 
inaccuracies are further increased due to the mentioned inaccuracies 


in the classification of the precipitation form. 


Air temperature 


Platinum resistance thermometers with a resolution of 0.1°C were 
mounted in radiation screens and the air temperatures were recorded 
every 30 minutes. The O°C level was tested by ice calibration 

twice a year, and the deviation seldom exceeded 0.1°C. Simonsen 
(1976) found fairly good protection against longwave loss as well 

as shortwave gain by using radiation screens of the type used in 
this investigation. Skartveit (1976) estimated the deviation between 
air temperature observed with an aspirated psyhcrometer and recorded 
by platinum thermometers in radiation screens of the same type 
during daytime in clear weather in snow-free conditions. He found 

a systematic deviation of 0.2°C and a standard deviation of 0.7°C. 
The low systematic deviation indicates fairly good performance 

of the radiation screens. Different sampling frequencies and the 
turbulent character of the air temperature are probably the main 
reasons for the standard deviation of 0.7°C. During snow covered 
ground conditions some higher deviations may occur, due to the 

high reflected radiation flux from the snow cover. 


The standard deviation of 0.7°C of the difference between the two 
types of temperature records referred to, indicates a standard 

error of 0.5°C for each of the record types. This estimate is adopted 
for the half-hourly records at the Dyrdalen-Frotveit area. For daily 
averages, this error obviously is reduced (to 0.1-0.2°C). 


Air humidity 


Relative humidity was recorded by Lambrecht Feuchtgeber 800 with 
an 11 cm long hair bundle placed in a radiation screen. The 
instruments were calibrated at the laboratory once a year, and 
routine checks against an Assman psychrometer in the field were 


a2) « 


performed. The accuracy of the corrected values are within t o¢ Ra. 
(Skartveit, 1976). If the standard error of the half-hourly values 
of the air temperature, T , is z 0.5°C, this will be approximately 
equivalent to a standard error in the vapour pressure, e , of 

+ 0.3 mb when T = 3°C . For daily means the error should be 


within z O.1 mb. 


Wind speed 


The wind speed was recorded by a 3 cup anemometer with a resolution 
of 115 m integrated wind travel (0.06 ms! for 4 hourly means), 

and a threshold speed about 0.3 ms! Skartveit (1976) tested the 
stability of the calibration factor and found it to be within + 14. 
The rotor is mounted between two horizontal solid plates, which 

may introduce spurious aerodynamical effects. During winter the 
rotor may seize due to snow and ice, but during snowmelt this 


should not be the case. 


At Dyrdalsvatn a 80 cm high fence surrounded the micrometeorological 
station to protect against grazing sheep and cattle. This may 
influence the wind speed at the 2 m level. The lysimeter was placed 
at a distance of some 100 m outside the fence. Thus, the wind field 
may be different at the lysimeter and the site of recording, but 

the local variation is probably small, due to the openness of the 
area on this scale (Fig. 2.3). The wind speed at the wooded site, 
Frotveit II, was recorded 2 m above the lysimeter. The wind speed 

in the wood was low, so the relative error is high, due to the 


threshold speed of the anemometer. 


According to the discussion above the errors of the wind speed 
values are difficult to quantify exactly. A rough guess is 
0.2-0.5 m/s for hourly values and 0.1-0.3 m/s for daily values. 


Radiation 


Global radiation and reflected solar radiation were measured with 
a Gorczynski solarimeter with Sonntag thermopiles. According to 
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Sonntag (1964), this instrument has an accuracy of z 3%. The 
thermopile can be considered independent of temperature, solar 
altitude, and azimuth (Radiation Yearbook, 1980). 


The instruments were calibrated each year, and the stability of 
the calibration factors was within 3%. An albedo of 60-70% and 
an inaccuracy of + 3% in each of the daily recorded incoming and 
reflected shortwave radiation terms produce a standard error 

of some 10% in absorbed shortwave radiation. The errors of the 
two fluxes are probably correlated, making this standard error 


even lower. 


Net radiation, Qn » 1s defined by, 
Qn = Q,(1-a) + Qy = Qn ° (3.143) 


where Qs is the global radiation, a the albedo of the surface, 
Qi y the incoming and QV 4 the outgoing longwave radiation. Assuming 
the emissitivity of snow, e4% 1 (Anderson, 1976, used the value 
0.99), then Q|4 = o18 » Where o is the Stefan Bolzmann constant 


(5.669-10 °Wm *k~*) and Th, ‘is the snow surface temperature. 


Qn was measured by Siemen Ersking net radiometers. Calibrations 

made at the Radiation Observatory in Bergen show that the factors 
given by the manufacturer should be lowered by some 10%. This 

result was also found by B. Holmgren, Uppsala (pers.comm.). At 
Dyrdalsvatn only the incoming total radiation was measured, 

by changing the lower lupolen hemisphere for an adaptor. The 
instrumental temperature, T. » Was also recorded, and wehave from (3.1. 


Q, = Qi tole - Qs (3.1.4) 


where Qn is the net radiation at a black surface of temperature, T. 


It should be reasonable to accept a 10% error in the average records 
of Qn - For daily values this error will be lower. | is very 
precisely determined, and Qs : = 3% . For average values of 


- 33 - 


2 4 


for oT.’ , and 150 wm* for Oe « 
1 9 S 


100 Wm’* for Qi) » 330 Wm 


then Q), = 280 wm7 2 
The accuracy of Qy when recording Qn by the more accurate 
Schulze net radiometer is z 2% (Paulsen, pers. comm.), which is 

6 wm 2 in our case. It should then be reasonable to put the 
accuracy of Qy (daily values), when using the Ersking instrument, 


to t 8 wm? (= 3%). 


with an error of z 10 Wm or less. 


The outgoing longwave radiation, oT) » is constant when the snow 
surface temperature, Ta » is at the melting point. When the air 
temperature, T, » < 273 K , a snow surface temperature 


parameterization given by 
T = 2+ T. - 273 (°K) (3.1.5) 


may be used (Harstveit, 1981). Since the snow surface is at melting 
point during most of the field measurements used, the error of the 
term oT, is small. The errors influencing the accuracy of Qn 
are then due to absorbed shortwave radiation (< 10%) and Qy 


2 and 280 Wm 
2 


For average values of these quantities (some 50 Wm “4 
a simple error analysis produces a standard error of = 10 Wm" 


for daily values of Qn 
Dew formation on the hemispheres may some times enlarge the errors. 


Extrapolated data of incoming longwave radiation, Qy » recorded 

in Bergen were used to fill gaps in the measurements at Dyrdalsvatn. 
Harstveit (1981) estimates an average decrease in Qy ef 24 wm” 2 
from the Radiation Observatory in Bergen to Dyrdalen based on the 
height gradient suggested by Kuzmin (1961). During 1980-82 Qy 

was recorded at Dyrdalsvatn. A comparison of the records from 

Bergen and Dyrdalsvatn during spring snowmelt (65 days) showed 


2 with a standard deviation of 11 Wm *, 


a mean difference of 17 Wm. 
Given the limited accuracy of longwave radiation measurements and 
the possibility of differing cloudiness between the two stations, 
2 mean difference is in reasonable agreement 


2 mean difference. The value 17 wm 2 is 


the measured 17 Wm 
with the estimated 24 Wm. 
used in this investigation. 
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3.2. The Energy Balance Model for Snowmelt. 


The complete energy exchange for a snowpack can be written: 


Qt Op = Wt Qy +O +O + O (3.2.1) 
where 

Qu = energy used in melting the snow 

Qy = changes in internal energy of the snow pack 

Qn = net radiation of the snow cover 

Quy = sensible heat flux from the atmosphere 

Qe = latent heat flux from the atmosphere 

Qo = heat flux from the ground 


Qe = energy gained from rainwater. 


In the case of snowmelt some of these terms are negligible. The 
snow pack may be composed of water in liquid and solid phase. For 

a snow pack which is not ripe or isothermal any energy supply will 
first be used to increase the internal energy of the snow pack, 

by raising the temperature or increasing the free water content. 
Analogously, energy loss will correspond to a decrease of the 
internal energy. During parts of the spring snowmelt season such 
variations occur througout the day. In this investigation only days 
where the snow pack can be considered ripe and isothermal, except 
for daily variations, are used. For daily totals Q, is therefore 


a negligible term. 


Heat is stored in the ground during the summer season. In winter 
this stored heat may be used in melting the snow. The heat flux, 
Qc » may be expressed by 

dT 


_ G =o 
Qe = Ke gq (Wm), (3.2.2) 
a alg 
where Ko is thermal conductivity and a temperature gradient 


in the surface layer of the soil. U.S.A.C.E. (1956) found this 
heat flux to be in the range of 0-4 wm, equivalent to O-1 mm 
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daily snowmelt. Granger (1977) found about the same values in 
Canadian prairie environments. This flux is accordingly considered 


negligible in our case. 


When rain of temperature Tp < 0°C falls on a melting snow cover 
the rain water will be cooled to 0°C and the released heat used 
in melting the snow. For a ripe and isothermal snow pack 


- ; “2 


where Pt is the rainfall rate (mm day!) and C., and Py are 


spesific heat and density of water, respectively. Thus, a daily 
total of 10 mm rain at a temperature of 5°C yields Qe = 2.4 wm? ‘ 
which would produce 0.6 mm melt water. Therefore, the term Qe is 


in general negligible. 


Accordingly, net radiation, Qn » andthe turbulent fluxes of sensible 
(Qu) and latent (Qe) heat from the atmosphere are the significant 
energy fluxes feeding the energy consumption, Qn » in snowmelt, viz. 


Wy = Wt gta » (3.2.4) 


which is the basic equation of our energy balance model. 
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3.3. Turbulent Fluxes of Heat. 


Here we will first show how the turbulent fluxes can be estimated 
from measurements of air temperature, air humidity, and wind speed. 
Net radiation will be discussed in a subsequent chapter. 


Sensible (Qu); and latent (Qe) heat is exchanged between the surface 
and the atmosphere due to vertical gradients in the air temperature 
and vapour pressure above the snow surface. Since almost 100% of 

this heat is transferred by turbulence, these are highly wind- 
dependent. During neutral stability the turbulent heat fluxes may 


be expressed by 


Qu = Pal ne pu Mtg! ; C3e3e1} 
Pally 
Qe 7—3 U (e,-e,) 3 (3.3.2) 


where Qu and Qe are defined positive downwards, Po is the 
density and c the spesific heat of the air, T is the air tempe- 
rature (°C), e@ is the vapour pressure (mb) and u is the wind 
speed at the level a(ms~!), Cy is the bulk transfer coefficient, 

which from theoretical considerations assuming idealized experimental 


conditions, turns out to be 


Z 
Cc, = «f/(In oA)? , (3.3.3) 
D ZL 
fe) 
where « is von Karman's constant (0.41), Z the instrumental 


a 
height (m), and Eo is the surface roughness (m). y is defined by 
.-* P 

p 


: a yal 
y¥ = 0.622 L, (mb K ) 5 (3.3.4) 


where Ps is the atmospheric pressure (mb) and L. is the latent 
heat ‘of sublimation (ak~'kg7!). 


The stratification above a melting snow cover is generally stable, 
and the equations (3.3.1) and (3.3.2) should be modified according 
to the Monin-Obukhov similarity theory, which is also based on 
idealized conditions. Following Anderson (1976), 
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Q, = f(u)(T,-T,) (3.3.5) 
Q- = sludle,-e,) » (3.3.6) 


where f(u) is the wind function defined by 


: ee 
flu) = Paty Cpl Ri. Rip) u (32.3.7) 


Ril. is the critical Richardson number, and Rip is the bulk 
Richardson number, defined by 

a o 
where the temperatures are given in °K, and g is the acceleration 
of gravity. The Richardson number reflects the ratio of the 
consumption of energy by the buoyancy forces to the rate of produc- 
tion by wind shear. The critical Richardson number is the value 
of Rip beyond which turbulent conditions no longer exist. 
Businger (1973) found this value to be 0.15-0.25. Brutsaert (1972) 
found that increasing evaporation and radiation losses increase 
Ri. . He concluded that there is no critical value, rather a range 
of values between 0.25 and 0.5, below which turbulence is very 
likely, and above which turbulence is improbable. 


When Ri. and Cy are known we may compute the heat fluxes 
Qu and Qe from data of the surface temperature, and the air 
temperature, humidity and wind speed at a level Z. . To determine 


C5 the surface roughness, Zo » must be known. 


During snowmelt periods, T) and e, are assumed to be O0°C, and 

6.11 mb, respectively, which should in general be valid except for 
some calm, cold nights. Usually Rip exceeds Ri. during such 
nights, and the theory fails. Due to fluctuations in the average 
wind field, turbulence will occur as soon as Rip falls below Ri 
This "bursting" process is described by Brutsaert (1972). The 
average turbulence during these nights, however, is small, and 


Qyt Qe 


wind speeds prevail 


is negligible 
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(Harstveit, 1981). For stations where low 


during the day the theory will not be suitable. 


The required horizontal homogeneity of the parameters and input 
data is not fulfilled most places in Western Norway due to the 


rough terrain. 


At Dyrdalsvatn the small scale variations (100-300 m 


scale) are small when the snow covers most of the ground, whicn 


was the usual situation during the periods of recorded snowmelt 


discussed here 


. On a 


larger scale, however, the terrain is rather 


rough (Fig. 2.2), and this may make the method less accurate. At 
Austlihnylla the scale of horizontal homogeneity is smaller (50- 
100m), which further increases the inaccuracy. A snow-drift 200 to 


300 m long and 
drift does not 


20 to 30 m wide usually builds up there, and this 
disappear until June or late May. At Frotveit the 


diverse vegetation makes the application of the theory rather 


doubtful. 


Accordingly, the question arises whether a simpler wind function 


may be substituted for (3.3.7). The following functions are tested: 
f(a) = a,(1-Ri,Ri “Ty ey (3.3.9) 
] ] Be “ ine 
F,(u) = au + bo , (Bate 10) 
ng 
and F,(u) = au : (3.3.11) 


Now, we recall 
rate, Qu > may 
fluxes 


are tnen optim 
(2(Q) 
wnere Qn 
and on 


is the measured and 
is the number of days. AS a 


that if net radiation, Qn » 1S measured, the melting 
be predicted from eq. (3.2.4) provided the turbulent 


(Qu ‘ Qe) are estimated from any of the wind functions 
(3.3+8)" (5.3. 11) 
and wind speea. 


and the recorded air temperature, air humidity, 


Tne aajustable constants of tnese wind functions 


ized by minimizing the residual error o, 


: ie ie ; (3.9.12) 


Quy the computed daily snowmelt, 


measure of the goodness of the 


a 


models, we may also use the efficiency of the model, R2 , defined by 


sp-2 (3.3.13) 
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where SD is the standard deviation of the recorded daily snowmelt, 


defined here by 
1 
- —_ 2 
Sh = t 2 (Qn Qn) } (3.3.14) 


with Qh as the mean of recorded daily snowmelt. 


The wind function Ff, (u) implies an explicit account of varying 
stability. The optimal value of a is related to the temporal 
average of surface roughness, zB Eq. (3.3.9) then corresponds 

to eq. (3.3.7), and, recalling @q. (3.3.3), it is seen that ay 

is given by 


2 
ay = aol = Pap —S— (am PK"!), (3.3.15) 
tin 4" 
0 
and from this we get, 
1 pe | 
«“fe¢ }°ea * 
2, © 2.6 =P (m) (3.3.16) 


In this way it is possible to calculate a temporal average of 
surface roughness from assumptions of Ri. » and records of air 
temperature, air humidity, and wind speed at a level Za net 


radiation and snowmelt rate. 


The simpler wind functions (3.3.10), and (3.3.11) may be inter- 

preted as weighted averages of actual stability conditions. The 

linear wind function (3.3.10) should be suitable also for wind- 

sheltered sites in that bo may be interpreted as an average of 
the turbulent activity for zero wind speed records. bo should 

of course be positive if adequately determined. 


3.4. Degree-Day Models for snowmelt. 


The degree-day (or temperature index) model in its most simple form 
is a linear regression model, stating that the daily snowmelt, Sin ; 
is related to the daily mean temperature, T , by 


m ae y 
T* (mm day 


Nn 
I 
= 
i 
f 
4H 
+ 
Il 
Vv 


) (3.4.1) 


oO 
1 
JA 


where the degree-day constant, k (mm day 'ec™!) » and the threshold 


temperature, T* (°C), should be determined by regression analysis. 
The model says nothing about the physical processes causing snow- 
melt, and is therefore more or less a blackbox model. The advantage 
of the model is its simplicity, and its success is due to the fact 
that the air temperature is more or less correlated to air humidity 
and net radiation, and hence to the terms Qn 5 Quy » and Qe in 

the energy balance model (3.2.4). This intercorrelation is different 
for different weather types, different seasons, and different terrains, 
however, and there is usually no correlation between wind speed and 
air temperature. The greatest disadvantage is that it has to be 
optimized for the specific area before it can be applied with any 
confidence. In areas exposed to wind and changing weather it fails. 
Ostrem (1974) clearly demonstrates this point. He found fairly 

good correlation between runoff and air temperature (r = 0.74 

for the ablation season 1971) at Austre Memurubre (61°33'N, 8°30'E), 
a valley glacier in the central part of Southern Norway sheltered 
against westerly weather. The Alfotbreen glacier (61°45'N, 5°40'E), 
covering exposed mountains near the coast of Western Norway, 
however, did not seem to respond significantly to variations in 

air temperature (r=0.05). The investigation of Hendrie and Price 
(1978) clearly shows that different degree-day constants have to 

be applied for deciduous and coniferous forests. 


Various improved forms of the degree-day method are in use today. 
The daily mean air temperature may be substituted by daily maximum 
or daily minimum temperature, and the model may be improved by 
taking physical principles into account. Otnes and Rastad (1979) 
suggest the equation 


T 


s, = [k,(1 - cos 200) 4 kolLt - TAT, (3.4.2) 


where I is the day number of the year. Such a seasonal dependence 
should appear due to variations in global radiation and albedo, 


ripening of the snowpack, and small scale advective effects during 
periods with partly snow-covered ground. Kuusisto (1978) suggests 
improvements by adding a radiation term to the model. 


Use of different degree-day constants during humid and dry weather 
may also improve the model. This is suggested by Anderson (1976) 
and Obled (1973). The snowmelt, Sm, may then be given by, 


k, (T-T*,) 9 T > TS 
» dry weather 


7 0 > | < T* 
il = aa. (3.4.2) 
ko (T-T*,) 5 Tt. ®& T*, 
» humid weather 
0 « T= TX, 


Where the cloud cover, or a precipitation indicator may be used 
in classifying the weather type. Roughly speaking, however, all 
the mentioned refinements of the simple degree-day model more or 
less implicitly draw in fragments of the energy balance model. 
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4, Results and Discussion. 


4.1. Observed Data. 


The deviation from the longterm average of the air temperature and 
the total of precipitation during the winter and the spring 1979-82 
are shown in Table 4.1. The most distinct deviations are that the 
winters (1979-80) were colder than normal, April (1979-81) drier, 
May 1979 colder and wetter, and May 1981 warmer and drier. 


Table 4.1. Deviation from the longterm average (1931-60) air 
temperature at Bergen-Florida and percent of the long- 
term average precipitation total at Bergen-Fredriksberg 


during winter and spring 1979-82. 


Winter season 


(Jan., Feb., March) | April May 
1979 “1.97 112 % 0%2°C 63 % =-205 % 196 % 
1980 ~0 , 68°C 63 % +0.4°C 57 % +0.4°C 74 % 
1981 -0.4°C 146 % =). 0"6 41 % +i2.5.0.7C 44 % 
1982 +0.6°C 120 % “0.2 °C 124 % -0.6°C 127 % 


In accordance with the low air temperature, very few snowmelt 
episodes occurred in Dyrdalen during the winters 1979 and 1980, and 
in 1979 the snow amounts were particularly large (Fig. 2.3). 


The periods of recorded snowmelt during which the snowpack could 

be classified as approximately ripe and isothermal except for 
diurnal variations are shown in Table 4.2. The spring snowmelt 
season in 1979 should be commented on. During most of April the 

snow was melting. From April 29 through May 6 snowpack again accumu- 
lated. This newly fallen snow ripened very rapidly, but the snowpack 
was strongly layered in the first half of May. During the subsequent 
days the whole snowpack quickly melted off. 


Diurnal averages of air temperature, humidity, wind speed, and 
totals of precipitation and snowmelt at Dyrdalsvatn in 1979-82 

are shown in Fig. 4.1. Also shown are the radiation fluxes and 
albedo of the snow surface. The mean melting rate was 11.4 mm day 
with a standard deviation of 15.3 mm day |. On May 21, 1979, there 
was ae i a melting of 110 mm(which requires an energy input of 
430 Wm “). 


] 


At Frotveit measurements were carried out within a deciduous wood 
and at an open site nearby. Fig. 4.2 shows the weather conditions 
during recorded snowmelt 1981 and 1982. Air temperature and humidity 
are nearly identical at the two sites. On the average the air 
temperature is 0.1° lower and the relative humidity 2% higher within 
the wood than outside. This promotes a corresponding increase of 
0.1 mb of the water vapour pressure. More significant deviations 
are found in the wind speed records, average values being 0.3 ms~ 
and 0.6 ms !, respectively, and by correlating daily values observed 
at the two sites, a correlation coefficient of only 0.58 appeared. 

A check against maps shows that for southerly winds (up-valley), 
which are often accompanied by precipitation, the wind speed is 
considerably higher at the open site. For down-valley winds, which 
often prevail during clear weather in early spring, the wind speed 
may be slightly higher at the wooded site. The open site is then 
sheltered by the wood. 


] 


For 12 days during March-April 1982 the average wind speed at Dyrdals- 
vatn was 1.4 ms |, and only 0.4 ms | at the woodless Frotveit site. 
The low wind speed values at Frotveit are due to the topography and 


vegetation conditions. 


At Austlihylla (Fig. 2.2, 2.4) air temperature, air humidity, and 
wind speed were recorded together with snowpack discharge (snow 
lysimeter). Precipitation data extrapolated from Dyrdalsvatn (pp. 27) 


were used. 
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Table 4.2. Periods of recorded snowmelt in the Dyrdalen and 
Frotveit areas during 1979-82. 


; ‘ Number 
Station Period of days 
Dyrdalsvatn April 1 to 30. 1979 30 

May Tt “89 21 1979 21 
Febr. 15 to 22 1980 
April 12 to May 1 1980 20 
May 1 to 9 1981 9 
Febr. 3 to 14 1982 12 
March 25 to April 7 1982 14 
April 15 to 22 1982 8 
122 
Austlihylla April 17, and 
April 23 to May 1 1980 10 
April 1 to 20 1981 20 
May 1 to 10 1981 10 
March 28 to April 7 1982 1] 
April 15 to 27 1982 13 
May 7 to 24 1982 18 
82 


Frotveit I (II) April 3 (4) to 
I : woodless site 14 (12) 1981 12 (9) 


II: wooded site Jan. 16 to 18 1982 3 
March 23 to April 5 1982 14 


- 45 - 


Table 4.3. Available (+) and missing (-) radiation records at the 
open (FI) and at the wooded (FII) sites at Frotveit, and 
at Dyrdalsvatn during the period of recorded snowmelt 
at Frotveit 1981. Qs denotes the incoming shortwave 
radiation, a the albedo, and Qn the net radiation. 


Station term 3/4 4/4 5/4 6/4 7/4 8/4 9/4 10/4 11/4 12/4 13/4 14/4 


Q.. + + + + + + + + + + + + 
FL aQ. + + + + + + + + + + + + 
Q + + + + it + = = - - - - 
N 
Q, - - - - = + + + + + + + 
FId aQ. - = = - 
Qy + + - + + fi + + + + + + 
Q, + + + + + + + + - - - - 
Dyrdals- aQ + + + + + + + + - « . 7 
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Fig. 4.1. Diurnal values of recorded data during snowmelt at 
Dyrdalsvatn 1979-82. 
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Fig. 4.2. Diurnal values of recorded snowmelt and meteorological 
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the radiation fluxes are estimated (Chapt. 4.2.1) due 
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Fig. 4.3. Diurnal values of recorded data at Austlihylla during 
| snowmelt 1980-82. 
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4.2. The Energy Balance of Recorded Snowmelt. 
4.2.1. Modelling Radiation Fluxes. 


The snow surface absorbs, reflects and emits short-wave as well 
as long-wave radiation. The net radiation of the snow surface is 


accordingly defined by 


Qn = Qo (1-a) + O47 Qi ry (4.2.1) 


where Qs is the global radiation, a is the albedo of the 
surface, Qry the incoming and Qs the outgoing long-wave 


radiation. 


Global radiation 


Global radiation, Qs » 1S an important and highly variable energy 
flux, varying regularly with season and hour, and more irregularly 
with cloudiness. The most precise determination of Qs is by 
direct measurements. However, estimates of snowmelt are frequently 
needed in areas where records of global radiation are missing. 


Global radiation data observed in Bergen, Qsp » and at Dyrdalsvatn, 
Qsp » during the periods of snowmelt discussed, are compared here. 
Daily mean values from 102 days yielded the regression 


Q Q 
T = 0.85(q°8) + 0.06, (4.2.2) 
ex ex 
where Q is the extraterrestrial global radiatian. The correla- 


ex 
tion coefficient is 0.94 and the mean values of (qe) are 0.40 


e 
for both sites. This indicates that the global radiation flux in 
this area is fairly invariant within distances of 10-20 km. 


Q Q 
The equation (4.2.2) predicts ae 0.70 when 7 = 0.75. The 
i Qex Qex 


average value of Ta, for all days with values exceeding 0.60 
D 


is 0.69, while oD is 0.63 when averaged over the same days. 


Qex 


Since the station at Dyrdalsvatn is more screened by the mountains 


this is to be expected. 


Q 
Eq. (4.2.2) also indicates that the highest values of L. during 
ex 


cloudy weather are found at Dyrdalsvatn. The average value of 


Q 

8 for all days with Qcp < 0.30 is 0.17, while the average 

a Qsp 

value of c= 0.19 for the same days. This can be explained by in- 
eX ; 


creased diffuse radiation due to multiple reflection between the 


snow cover and the clouds. 


In accordance with these systematic differences, the standard devia- 
tion of the data from Bergen is 0.22, and 0.20 for Dyrdalsvatn. 

Some of the deviation from the line 1:1 by eq. (4.2.2) may, however, 
also be due to differences in cloud cover at the two sites. 


The cloud cover and the relative sunshine duration may both be 
used for estimating global radiation. These two variables are 
closely related. For the 102 days used for determination of eq. 
(4.2.2), the regression 


= |.18 (1-6) - 0,02 » YF = 6.91 (4.2.3) 


is found. C is the cloud cover, defined in this chapter as a 
fraction of unity, and z is relative sunshine duration, defined 

as the ratio of actual to possible sunshine duration (all data from 
Bergen-Florida). Some of the unexplained variance may be due to 
inaccuracies in observed cloud cover. The slope of the line (4.2.3) 
deviates from unity, which is probably due to systematic differences 
during fine weather. Thin clouds, or clouds not covering the sun 
will lower the values of 1-C , but not those of < . Besides, C 

is reported as 1/8 at the slightest trace of clouds, partly explain- 
ing why (4.2.3) predicts é ~ 1.0 at C= 1/8 

Different authors (Kimball (1914), Prescott (1940), Penmann (1948), 
McCay (1970), Bayne (1978)) found good estimates of relative 

global radiation, ae from é or C . Skartveit (1976) used the 


; ex 
regression model 


Qs 


ex 


S 


where C. = é » and S} > So >» and S3 are regression coefficients. 
This model is optimized on a daily basis for the 102 days used for 
determination of eq. (4.2.2), using both ‘ » and 1-C for C. 
Results from the analysis made by Skartveit (1976) and the coeffici- 
ents found in this investigation are listed in Table 4.4. Note 

that s/S and C€ are observed in Bergen only, while the computa- 

tion is made for global radiation recorded both in Bergen and at 
Dyrdalsvatn. From the table it can be deduced that me is pre- 
dicted higher by using the present coefficients than by using those 
found by Skartveit (1976). The difference is, however, hardly statis- 
tically significant. When using global radiation data from Dyrdals- 
vatn, a slightly lower correlation is found, as expected due to the 
distance between the sites. The regression coefficients differ 
somewhat, in accordance with the differences in global radiation 
(4.2.2). Use of relative clearness, 1-C , for C. yields a some- 
what lower correlation. This is partly due to the mentioned inaccura- 
cies in the cloud cover data, and partly due to z being a more 


adequate physical indicator. 


We may conclude that rather good Qo-estimates can be derived from 
cloud cover or relative sunshine data. This is a fact of practical 
importance since the cloud cover is observed at weather stations 
on a routine basis, so that fairly reliable data are available 

for most areas. The equation (4.2.4d) (Table 4.4) is then recommended. 


For comparing global radiation at Dyrdalsvatn and Frotveit (open site) only 

8 days are available. The mean values were Qsp = 119 Wm™é , and Qc = 116 Wm 
with a correlation coefficient» r = 0.96. By excluding a single day with fog 
at Frotveit (no fog at Dyrdalsvatn on that day), Qsp 


2 


112 Wm * , 

Qs] = 114 wim” 2 , and r = 0.99. This confirms the high quality of our data 
and shows that global radiation records from Dyrdalsvatn can be used as high 
quality estimates for the open site at Frotveit, FI , when data are missing. 
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Snow albedo 


The albedo of snow is generally high, but is shows quite substan- 
tial and very important variations. The crystal structure, the 
liquid water content of the snow cover, and litter, dust, etc. 

at the surface all influence the albedo. These parameters change 
with time, and Fig. 4.4 (from U.S.A.C.E., 1956) exemplifies how 
the resulting time dependency of snow albedo may be parameter- 
ized in terms of time since last snowfall Or in terms of accumu- 
lated daily maximum air temperature. A regression analysis between 
the albedo, a , and the time since last snowfall (days) during 
spring snowmelt at Dyrdalsvatn 1979-82 (91 days, days with snowfall 
excluded) yields a correlation coefficient, r , of only -0.24. 

We now define the time, t » Of open air exposure of the snow at 
the surface by 


n+ 1 - 3T. when the old surface is 
exposed to open air 
E = (days) (4.2.5) 


n. + 1 when the old surface is 

covered by a new snowfall, 

where n is the number of days since the last day of snowfall 
before spring snowmelt, and n. the number of days since the last 
snowfall. T. is the number of days when a new snowfall is cover- 
ing the old surface, estimated from totals of snowfall and a first 
order estimate of snowmelt. The definition is illustrated in 

Fig. 4.5. It should be noted that this is a very rough parameter- 
ization of the surface conditions. Nevertheless, the correlation 
coefficient between a and t was -0.50 for the 91 days at 
Dyrdalsvatn 1979-82, indicating that the sudden changes in the 
albedo due to the covering and uncovering of the old surface 
affects the albedo Significantly. 


The relation between a and t is illustrated in Fig. 4.6. The 
figure indicates a more rapid albedo decrease for the first days 
after a snowfall. By exchanging t with Int =. the correlation 
coefficient, r , increases to -0.55. As can also be seen, the 
albedo decrease with time is slower than Suggested by U.S.A.C.E 
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VARIATION IN ALBEDO WITH TIME 
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Fig. 4.4. Variation of albedo (after U.S.A.C.E., 1956). 
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(Fig. 4.4). However, that result was based on data from Centra 
Sierra, situated at about 39°N, and with a continental climate. 
Due to the lower latitude, the crystal transformations are 


expected to work faster in that area. 
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Fig. 4.5. Illustration of the definition (4.2.5) of the time, 
t (days), of open air exposure of the snow at the 


surface. 


Fig. 4.7 shows the regression line between the albedo and the 
relative clearness, 1-C . The correlation coefficient is -0.42 
which is rather high considering the neglection of variations 
in snow surface age. An increase in snow albedo with increasing 
cloud cover is also reported by Kondratjev (1969, pp. 426), 


56 


"28-6261 uzeAS|eputg 
ye zpawmous BuLuds Bbulunp * » Sopaqyre ayy pue (€G *dd) ade suns 
au 2e@ MOUS B3Y2 JO AauNsSodxsd ule Uado JO * Y S9WLZ DYY UdSM}9q UOLZeLAY “O°y “BLY 


RO SS Bio gone 


(SS°0- =4) 080+4u) 400 


I 


57 


*uabueg ye paruasqgo * J-| ‘Sssauueald aALyelau pue °Z28-6/61 
‘uzeastepukg ye ylauwmous Bbuluds Bulunp * » fopaq,e ayy uaamizaq UOLZELOY “/*p “BLY 


L0 


80 


- 58 - 


"(9G6L S°F°I WSN 4aqse) 
Mous BuLlqLaw 40 AZLALZIALJau Leuqoads arzewLxouddy 


SNOYIJIW “HLONST 3AVM 
9°L TI Zh Ol 8°0 90 70 20 


(SE6L) XNDABQ Ja4sy xX 
(L76L) JoulDs, Jasjv O 


“9° 


v 


“Bid 


LNI3YId * ( O0g391V) ALIAILIII4 34 


« 59 = 


and Liljequist (1956-57). The line in Fig.4.7 should perhaps be 
curved towards high values of a for dense cloud cover, a 
tendency also reported by Kondratjev (1969, pp. 426). The scatter- 
ing process in clouds is only slightly spectrally dependent in the 
shortwave band (Robinson, 1966, pp. 173), and the transmission 
through the pyranometer glass bulb is nearly independent of the 
wavelength between 0.3» and 2p (Robinson, 1966, pp. 265; 
Sonntag, 1964). 


The spectral albedo of snow, however, increases with decreasing 
wavelength (Fig. 4.8). Multiple reflections between the snow 
surface and the clouds therefore simultaneously enhance the global 
radiation and cause a "blue-shift" in its spectral distribution. 
Such a "blue-shift" in turn enhances the integral snow albedo, 
given the above-mentioned correlation between spectral albedo 


and wavelength. 


A multiple correlation analysis albedo and the relative clearness, 
1-C , and In t (pp.53) yields a regression model (91 days) for 
the albedo, ao (in fraction of unity): 


a= -0,13(1-€) - 0.05 In € + 0.87 , (4.2.6) 


where C is observed in Bergen. The correlation coefficient is 

R = 0.75 , which means that 56% of the variation of a is explained. 
A slightly better model is found by replacing relative clearness 

with relative global radiation, = » where Msp is global radia- 


tion observed at Dyrdalsvatn: 


; Qsp 
a = -0.18 ~— - 0.05 In t + 0.90 , (4.2.7) 
Vex 

where R=0.77. Both models predict a 12% decrease of the albedo 
10 days after last snowfall if older surfaces are not exposed, 
and a corresponding decrease during the transition from overcast 
to clear weather. The unexplained variance is rather high, as should 
be expected considering the many factors which influence the albedo. 
If possible, these models should be optimized when used at other 


sites. This is especially important when regional climate, height 
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above sea level, snowmelt season, vegetation or pollution of the 
snow surface deviate significantly from those of Dyrdalen. 


By comparing the albedo recorded at Dyrdalsvatn and the open site, 
FI , at Frotveit, the values at FI are found to be some 10% 
lower than at Dyrdalsvatn, on the average. The correlation coeffi- 
cient is 0.79, but the sample is not large enough to allow a 
discussion of the difference. Falls from the vegetation in the 
area, and the lower altitude (higher temperature) probably account 
for most of the albedo decrease. When missing at Frotveit, Or 

is therefore extrapolated from Dyrdalsvatn by subtracting 10%. 

The exact value of the albedo at the wooded site, ArT» is some- 
what less critical due to the lower incoming short wave radiation, 
and was not recorded. This value is put at 10% lower than Or > 
since the surface within the wood appeared darker to the eye due 


to more fall from the trees. 


Longwave radiation 


In contrast to the case of shortwave radiation, the absorbitivity 

of a snow cover is approximately unity for longwave radiation. 

This fact enhances the importance of downward longwave radiation, Q , : 
relative to that of global radiation in connection with snowmelt. 
When not recorded, incoming longwave radiation during cloudy 

weather is frequently estimated from formulae of the following type: 
ey (4.2.8) 


1¢Li Ce Lge 


2 


2 


GW, = Gry | 
where C is cloud cover (Angstrom (1913), Bolz (1949), Paulsen 
and Torheim (1964), and others). The coefficients L,' and Ly! 
are sometimes assumed to be functions of cloud type (Monteith, 1975, 
op. 37): QoLy is the atmospheric radiation during clear weather 
at the same air temperature (and humidity). For Bergen, Paulsen 
and Torheim (1964) arrived at an equation for recorded hourly 
values during clear weather, given by Paulsen (1967) as 


= 0.895 o T’ - 29.3. (Wm 2), (4.2.9) 


QoLe 


where T is the air temperature (°K) , and o is Stephan- 
Bolzmann's constant. The correlation coefficient was r = 0.949 
Now, the formula (4.2.8) prescribes Qt 4 to be multiplied by 

a cloud factor to obtain actual atmospheric radiation, Qy 
Paltridge and Platt (1976, pp. 138-140) argue that the radiation 
outside the atmospheric window (8-13um) is of a black body 
character, and is independent of the cloud cover. The extra radia- 
tion from clouds must therefore occur in the window region. They 


then suggest the equation 


+ (1-0.7)e, T C (ls (4.2.10) 


Wy = Bory 
where Oot + is the radiation outside the window region, Eo is 
cloud emissivity (0.3 for cirrus; 1.0 for low clouds), and 
Tz cloud base temperature. The factor (1-0.7) is the approximate 
part of blackbody radiation which occurs between 8 and 13 um. 
Fig. 4.9 shows an approximately linear increase of the incoming 
longwave radiation, Qiy with the cloud cover, C . A regression 


model like 


2 ieet” + & 


1 (4.2.11) 


Qy 
should then be reasonable. at is blackbody radiation at air 
temperature, and C is a cloud variable which may include cloud 
type and height of cloud base. Ly : Lo » and L4 are the regress- 
ion coefficients. 


Table 4.5 shows wie eon) es of the multiple regression analysis 
using C, ¢* » and (qa) for Cy. C, 
150 monthly mean values (1965-1979) from Bergen. When the 


coefficients of the various formulae is compared it should be 
2 
) 


= C is also applied to 


remembered that when C increases from O to 1 , (sf, 
declines from about 0.5 to O . The predicted contribution to 
atmospheric radiation, Qy » from clouds during completely over- 
cast weather should then be 80, 66, 59, and 71 (wm?) for the 

four versions, respectively. This is reasonably close to the 

value 60 Wm 2 found by Paitridge and Platt (1976, pp. 140) 

for Aspendale, Australia, " a fairly typical mid latitude station." 
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Fig. 4.9. Relation between cloud cover, C , from Bergen, and 
incoming longwave radiation, Qy » at Dyrdalsvatn 
during 65 days of spring snowmelt 1980-82. 
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In Fig. 4.10, predicted Qy values for Ml, M2, M3 and M4 during 
clear weather (C=0) are plotted against the well known formula 
derived by Swinbank (1963): 


= 8.31 * 10° (wm?) (4.2.12) 


Qy 
and the formula found by Paulsen (1967) (M5). The formula based 
on the large set of monthly mean values from Bergen predict Qi y 
values under clear skies very close to those found by using the 
Swinbank formula, and notably lower values than by using the 
Paulsen formula. All the three model versions MI, M2, and M3 also 
predict lower clear weather values of Qry than found by using 


the Paulsen formula. 


The somewhat low increase in Qry with increasing air temperature 
found in Ml, M2, and M3 may be due to the relative narrow tempera- 
ture interval used (276.5 * 3.0°K) compared to the rather high 
variance in the cloud cover data (Fig. 4.9). Therefore, these 
equations should not be applied too far outside the temperature 
interval mentioned. In that respect M4 should be more reliable. 

It is highly noteworthy that this equation yields a cloud cover 
term in accordance with Paltridge and Platt (1976) and a clear sky 
term in accordance with Swinbank (1963), a fact that strengthens 


our confidence in that particular model. 
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Qy):5.31: 10° T§( Wm?) ( Swinbank, 1963 ) 


: 0.70 oT*+ 6.2 (Wm?) 
: 0.75 o T4+ 7.5 (Wm?) 
: 0.73 oT*+ 18.9 (Wm2) 
: 1.02 ¢T4- 91.5 (Wm-2) 
: 0.895 o T4 - 29.3 (Wm) 


190 200 210 220 230 240 250 260 270 280 290 300 310 
Qi y: Mi (Wm72) 


Fig. 4.10. Estimated longwave radiation, Qi y for clear weather 
using the model versions Ml to M4 (see table 4.5), 
and the Paulsen formula (M5) plotted against the 
Swinbank formula. 


The effect of a leafless, deciduous wood canopy on downward 


radiation. 


Any wood canopy affects all energy fluxes involved in snow melting. 
The observations at Frotveit make possible a rough quantification 
of how, in particular, a deciduous wood modifies the input of 


radiant energy to the melting process. 


The downward. longwave radiation, Qry , over a "black" surface of 


temperature, T, » May be expressed from eq. (4.2.1) as: 


= Qy - (1-a)Q, + oT 4 (4.2.13) 


QLy (e) 


where Qn is the net radiation and Qc is the global radiation 
at the snow surface, while a is the albedo and o is the 


Stephan-Bolzmann's constant. 
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Fig. 4.11. Recorded net radiation, Qn » at Frotveit (woodless (I), 
and wooded (II) site), and incoming longwave radiation, 
Qy » at Bergen-Florida (B) and Dyrdalsvatn (D) during 
the clear night of April 6-7 1981. 
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At Frotveit Qn ‘ Qs and a were measured at the woodless site, 
and Qn and Qs at the wooded site. Due to fall from the trees 
the albedo is lower within the wood than outside. As discussed 
earlier, this difference is assumed to be 10%. Moreover, T=273°K 
whenever the snow is melting, and Qy is readily obtained from 
eq. (4.2.13) insuch cases. A comparison of daily means thus derived 
for the woodless Frotveit site and corresponding values measured 

at Dyrdalen, yields a correlation coefficient of 0.97 and a mean 
difference of only 0.9 wm %, Considering the different ways by 
which these Qry are obtained at the two sites, this data comparison 
is in remarkably good agreement with the small intersite difference 
expected (< 5 wm” 2) from the small horizontal and vertical distance 
between the sites. In the case of missing data at the woodless 
Frotveit site, Qy values measured at Dyrdalsvatn are therefore 
used without corrections. When snowmelt does not occur, however, 
the surface temperature Ts may well drop substantially below 273°K. 
The cloudless night between 6 and 7 April. 1981 provides one, not 
extreme, example of this. The available snowmelt data indicate 

that the surface temperature stayed at 273°K until 20 hr CET in the 
evening of 6 April at both Frotveit sites. Wy at that hour is 
therefore known at both sites from eq. (4.2.13). Moreover, the 
recorded Qry follow quite identical patterns at Dyrdalsvatn and 
at Bergen during this night (Fig. 4.11). This observation strongly 
corroborates that the changes in Qy during the night are the same 
at the woodless Frotveit site as Bergen and Dyrdalsvatn. Further- 
more, the excess Qi y at the wooded Frotveit site (relative to Qi y 
at the woodless site), due mainly to the downward "window" radia- 
tion (8-13 um wavelength) from the tree canopies, is assumed to 
remain at its 20 hr CET value (42 Wm” 2) throughout the night. This 
last assumption is reasonable since the downfacing canopy parts 

are close to the air temperature, which varies only slightly 


during the night (Fig. 4.12). 


By the above assumptions we have estimates of Qi y throughout the 
night at both Frotveit sites (Fig. 4.12), and the corresponding 
surface temperature, Ts > 1S readily obtained from eq. (4.2.13). 
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4.12. Recorded air temperature, Li ,» and estimated incoming 
longwave radiation,.Q. | » and snow surface temperature, T 
at Frotveit during the clear night of April 6-7 1981. 


O 


I : Woodless site 
Il: Wooded site 


It is seen (Fig. 4.12) that T, at the woodless site drops rapidly 
before midnight and then levels off at about -4°C, while the air 
temperature 1,8 m above the snow surface remains at about -0.5°%C 
throughout the night. The existence of such a surface inversion is 
qualitatively corroborated by the observation that the air tempera- 
ture at Austlihylla (Fig. 2.2) this night was above 0°C while it 
dropped to -7°C at Dyrdalsvatn. At the wooded site, on the other 
hand, the estimated surface temperature dropped only slightly 

below 0°C during this night, well in accordance with the common 
observation that during cold spring nights the snow surface crust 
forms far more rapidly in woodless areas than within a wood. 


The difference in downward radiation at the wooded and at the wood- 
less sites may be interpreted in terms of the canopy porousity 
(Hendrie and Price, 1978) or the sky obscuration factor (Perthu, 
1982). The canopy porousity, Pe » for isotropic downward shortwave 
radiation may be measured accurately during overcast conditions. 
During two overcast days the ratio Pe between the mean global 
radiation (Qo) at the wooded site (II) and the woodless site (I) was: 


Q 
p: - SI] _ 9 5 (4.2.14) 


Assuming again that the downfacing canopy parts radiate at air 
temperature, T. » and that the downward longwave sky radiation is 
isotropic, we derive an alternative measure Pa of canopy porousity 
from the following equation, relating Qiy at the wooded (II) and 
the woodless (I) sites: 


(TPot eT a + Po'Q yr = Qyry ada 
During the cloudless night 6 - 7 April 1981 this Pe varied slightly 
about the value 0.47. Since the difference between black-body radi- 
ance at air temperature (“canopy radiance) and clear sky radiance 
increases with decreasing zenith angle (Kondratjev, 1969) the above 
Re is significantly biased towards the canopy porousity just over- 


head, which may by accident deviate considerably from an average value. 


- 70 - 


Since the sky radiance is far more isotropic in cloudy weather the 
Following mean value is chosen, as a compromise, for the average 
porousity for downward longwave radiation at the wooded site: 


tele a 3tP Ped = 0.53 ‘ (4.2.16) 


It should be mentioned here that this canopy porousity is a single 
spot value of an obviously variable quantity. 


A substitute of a canopy porousity of 0.53 into eq. (4.2.15) yields 


4 


= 0.4707. + 0.53 Qa] : (4.2.17) 


hart 
If not available, Qa] may be estimated from the empirical formula 
found for Dyrdalen (Table-4.5): 


4 


= 0.7307. 118 (~~ 


Cyr = ep 
which yields by substitution into eq. (4.2.17): 


Q _ 
4 62 (99)? 5 ay (wm?) © (4.2.19) 


= 0.860T 
. Qax 


Wert 
By comparing tnese formulae for Qry at the wooded (I1) and the 
woodless (I) sites, it is seen that the tree canopy reduces the 
sensivity of Qiy to variations in real cl oats (parameterized 
in terms of the relative global radiation Tex ) and enhances its 
sensivity to air temperature. This merely reflects the analogy 
between the tree canopy and a constant partial cloud cover at air 


temperature. 


The period of simultaneous registrations of shortwave radiation at 
I and If (Table 4.3) includes cloudy days as well as clear days. 


The line of regression is given by, 
oj 


Qst] = 0.71 Qc] * 4.7 (Wm 3 (4.2.20) 


and the correlation coefficient is 0.994. The regression coefficients, 
reflecting the shading from the wood are, however, spot values 


valid only for the actual site. The average yalues of Qo] and 
Qo7z are 99 wm” and 66 Wm < 
reduces the downward solar radiation and simultaneously enhances 
the downward longwave radiation. During 6 days of parallel record- 
ings, the net radiation, Qn , was on average 1/7 Wm 2 higher at the 
wooded site than at the woodless site, which is qualitatively in 
agreement with the observations reported by Petzold and Wilson 
(1974) and by Hendrie and Price (1978). This difference between the 
wooded (II) and the woodless (I) sites, however, varies signifi- 
cantly with cloudiness. The following regression equations relating 


» respectively. The tree canopy thus 


daily mean values were found: 


\Qy = QNTI ~ Quy = 214 1-C) + 5 (wm *) (4.2221) 
AQy = Q - Qyy = 36(05)2 # 5 (wm?) V4.2 e228) 
N > XNIT NI To, 2. 


with correlation coefficients of 0.67 and 0.79, respectively. 

Here C is fractional cloud cover observed in Bergen, while Qs 

is the global radiation measured at Frotveit or Dyrdalsvatn. The 
lower correlation coefficient in the former case is thus at least 
partly explained from the fact that the cloud cover observed in 
Bergen may occansionally deviate from Frotveit/Dyrdalen. These 

two regression equations are based on 26 days, during more than 
half of which radiation records are missing at Frotveit = and 
therefore extrapolated from data at Dyrdalsvatn as described above. 


According to eas. (4.2.21 - 4.2.22) net radiation is higher within 
than outside the wood, and this difference is most pronounced in 
cloudless weather. The contributions to this pattern from the various 
radiation fluxes are illustrated in Fig. 4.2 and Table 4.6. It is 
seen that the lower albedo within the wood counteracts the shading 
of global radiation. The increase in net radiation within the wood 
is thus largely controlled by the enhancement of downward longwave 
radiation due to "window" radiation from the tree canopies. 
Consequently, this net radiation increase is enhanced by reduced 

sky radiance within the "window" region, i.e. by reduced cloud cover 
in accordance with eqs. (4.2.21 - 4.2.22). 


Table 4.6. Global radiation (Qo); surface albedo (a), net shortwave 
radiation (K*), blackbody radiation (or, "1s downward 
longwave radiation (Q) ,)> and net radiation difference 
(Qu) between wooded (II) and woodless (I) site at 
Frotveit. All data refer to close to snow surface level. 

: I II 4 I Il 
Periag a Sr ST SF, 7, xa Sia Ben 
6 days 
Qn recorded 136 92 0.59 0.49 56 46 315 276 303 17 
26 days 
Qn partly 86 57 0.62 0.52 33 28 315 302 317 11 


modelled 


4.2.2. Turbulent Fluxes of Heat, and Model Fit to observed Data. 


The energy consumed in snowmelt, Qn » may now be estimated from the 


equation (3.2.4), i.e. 
Qu = Qy + Qy + Qe; 


where the net radiation, Qn » 18 recorded or estimated as described 
in Chapt. (4.2.1), and the turbulent fluxes of sensible (Qu) and 
latent (Qe) heat are estimated from the equations (3.3.5) and (3.3.6). 
Half-hourly values of wind speed, u , air temperature, Ty » and 
Vapour pressure, e, , recorded at the height, Z. » above the snow 
surface (1.3 m on the average) are used in calculating daily totals 
of the heat fluxes. The three wind functions (3.3.9), (3.3.10), and 
(3.3.11) are applied, and the optimal values of the adjustable 
constants (a, b) in those functions are found. This is done by com- 
paring daily totals of computed and recorded snowmelt and determining 
the set of (a, b) values which yields the lowest residual error, 

On (3.3.12), 

Results from Dyrdalsvatn snowmelt station. 


The optimal values of a and b and the corresponding model 
efficiency, R2 (3.3.13), are given in Table 4.7. The three wind 
functions yield nearly identical efficiencies. Therefore, only the 
wind function, fF, (u) » (3.3.9), which contains a stability correction 
factor (1-Rigki, |), and the more simple and widely used linear 
function, f,(u) » (3.3.10) will be further discussed. 


Following Anderson (1976), Ri, = 0.4 is applied in Ff, (u) » This 


yields by optimalization a, = 4.05 gm7 3x7! , and, from eq. 


(343516) Z,.= 0.08 cm . Optimizing with ein = 0.2 gives 

a. © 4.40 am °K" | and he = 0.11 cm , indicating that the model is 
not very sensitive to Riv. Sverdrup (1936) arrived at Za = 0.23 cm by 
measuring wind profiles above a melting snow cover at West Spits- 
bergen, while Liljequist (1956-57) found values between 0.01 and 


0.1 cm above an Antarctic snow field. Kuzmin (1961) suggests 
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Z,=0.05 cm for a stable snow cover deeper than 10-20 cm. The 
investigation of Anderson (1976) indicated varying but evidently 
decreasing values of Zo through the melt season, with a mean 
value of 0.15 cm. Accordingly, a value of 0.08 cm seems to be 
reasonable for Dyrdalsvatn during snowmelt. 


Table 4.7. Model fit and optimal values of the adjustable constants 
(a, b) in the three wind functions during snowmelt at 
Dyrdalsvatn 1979-82. 


Residual Model Standard 
error efficiency a. Ds deviation 
Wind Function ov R2 1 
mm day! win? gm 3K"! mm day !Wm 2 
f,(u)=a, (1-Ri "Rig)?u 6.94 26.7 0.793 4.05 7 
fo (u)=a,utb, 6.88 26.5 0.797 3.65 0.0(Wm 7K!) 
F3(u)=azu'3 6.94 26.7 0.793 4.05 0.95 
Recorded snowmelt 15.3 58.8 


The standard deviations in Ty » U, and Z. lead to a standard 
deviation of 10% in the stability correction factor (1-RigRi.”'). 
From Table 4.7, however, it can be seen that the introduction of 
this variable factor does not improve the model. Possible improve- 
ment may be masked by variation in Zoe and/or noise in the records 
of net radiation and snowmelt. The deviation from ideal site 
conditions may also strongly affect the stability factor, which is 


explicitly related to ideal conditions. 


The optimal wind functions adopted for the Dyrdalsvatn snowmelt 


station are: 


fy(u) = 4.05(1-RigRi,”!)2u (wm72) (4.2.23) 
where RT = 0.4 > and 
fo(u) = 3.65 u, (wm™*) (4.2.24) 
where the wind speed, u , is given in ms”! - During neutral conditions 


Rip = 0 , which yields Ff, (u) = 4.05 u 
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Since (4.2.24) is developed from all observations at Dyrdalsvatn, 
this indicates that stable conditions above the melting snow cover 
Suppress the average turbulent exchange by some 10% from the 
expected value during neutral conditions. 


The topography of the model efficiency, R2 , when f,(u) is applied, 
is illustrated in Fig. 4.13. The optimal area is rather wide, 
indicating that many combinations of (a, b) values are suitable. 

By extrapolating the "ridge" to a=0 , the wind function degenerates 
to a constant. The efficiency, R2 , however, then is notably lower 
than the optimal value. It therefore can be established that varia- 


tions in the wind speed is important to model. 
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Fig. 4.13. Part of the R2 topography for the wind function 
Ff, (u) = au +b for data from Dyrdalsvatn during 
snowmelt 1979-82. 


For the linear wind function, Ff, (u) = au +b , Harstveit (1981) 
found the following (a, b) values: for April 1979 (3.1, 0.0), 
for May 1979 (3.85, 1.9), and for April 1980 (3.65, 1.15). 


These values are slightly biased towards higher values, due to a 
systematic error in the incoming longwave radiation values, which 
has now been adjusted (Chapt. 3.1). Though there is some variation 
from one period to another, the constants are comparatively stable. 


The figures 4.14 and 4.16 illustrate the fit of the energy balance :: 
model to recorded snowmelt when the optimal linear wind function 
(4.2.24) is applied. The fit is generally good. During the days 
May-12 .=°15,. 1979, however, large deviations are found. In 
the first part of this month, the old, coarse-grained snow pack 
was covered by a finer-grained layer, and the snow pack thus was 
strongly layered. Some of the melt water produced during the days 
May 1 - 11 was probably stored in the snow pack since the vertical 
flow of melt water is impeded by such a layering (Jordan, 1983). 
During the days of rapid snowmelt, May 12 - 15, the stored melt 
water percolated together with actual rain and melt water, thereby 
giving too high recorded melt rates. By excluding those days, the 
efficiency, R2, increases from 0.80 to 0.84. 


Results from Frotveit. 


The topography and local variations in the surface roughness at 
Frotveit make the application of the theory for deducing the wind 
functions rather doubtful. Nevertheless, an attempt was made to 
compute the constants in the linear wind function from the data set. 


The optimal function was: 
(u) = 3.1 u + 2.3 (wm 2K7!) (4.2.25) 
for the woodless site, I, and 


(u) = 5.4u 4+ 1.9 (um 2K7 4) (4.2.26) 
for the wooded site, II. The wind speed, u, is given in a”, 

The efficiency, R2, is 0.71 at the woodless site and 0.83 at the 
wooded site. The difference could partly be explained by the dif- 
ferent wind speed patterns at the two sites: At the wooded site 
the wind speed remained very low through all the days of snowmelt, 


while at the woodless site it was relatively high during humid 
weather conditions (Table 4.10). This could result in a more 
variable turbulence mechanism at the open site. In addition, the 
net radiation is greater within the wood due to longwave radiation 
from the trees (Chapt. 4.2.1). The radiation estimates may be more 
precisely determined than the turbulent heat fluxes, and thus the 
total energy balance estimates may be of higher quality within 

the wood than outside. The samples are small, however, and the 
difference in R2 between the two sites may be insignificant. 


The optimal wind functions, Ff, (u) and Fry (u) » are used in the 
further calculations. At a first glance, these functions may seem 
rather different. Both indicate, however, that the turbulent 
activity does not vanish for the very low wind speeds which pre- 
vailed at Frotveit during the recorded snowmelt periods. Businger 
(1973) suggests a "bursting" process during very stable conditions. 
The air near the snow surface is cooled, and small wind speed 
fluctuations mix the air when the layering becomes too strong. 


The function fr, (u) deviates from fF, (u) mostly through the 
higher value of the slope constant, a. This could be explained 

by higher surface roughness values. However, the model fit is rather 
insensitive to variations of this adjustable constant due to the 
very low wind speed prevailing in the wood. By exchanging (4.2.26) 
for (4.2.25) at the wooded site, R2 drops only from 0.83 to 0.79. 


Results from Austlihylla. 


For Austlihylla, the optimal linear function is 
_ =P, = | 
flu) = 3.1 u (Wm °K ~) (4.2.27) 


where u is the wind speed (ms~!), This function yields somewhat 
lower values of the heat fluxes than the optimal function at 
Dyrdalsvatn during identical weather conditions. The model effici- 
ency, R2 = 0.47 , which is a markedly poorer: model fit than was 
found at Dyrdalsvatn. 


Due to the large amount of snow at Austlyhylla a time lag between 
actual and recorded snowmelt/rain amount should be taken into 
account also on daily basis. The model improvement when using an 
optimal time lag (4 hours) between actual snowmelt/rain input and 
recorded discharge is small (Table 4.8), however, and cannot explain 


the low model efficiency found. 


The data set is also separated in two parts according to the amount 
of daily precipitation, P : P < 5 mm (category I, 60 days), and 
p > 5 mm (category II, 22 days). 


Table 4.8. Model efficiencies, standard deviations, and mean values 
of modelled and recorded snowmelt at Austlihylla during 


1979-82. 
Model Standard Mean 
No Classi- Efficiency Deviations Values 
of days fication Model R2 mm day ‘Wm 2 mm day ‘wm? 
82 os Energy Balance 9 a7 11.6 44.7. 13.3 51.2 
FluJe3.1 u . , , , , 
Energy Balance 
Tlij=3.1 wu 
82 — and & heuine 0.50 11.5 44.3 13.3 51.2 
time lag 
60 P<S mm Energy Balance 4 ¢5 12.1 46.6 12.8 49.3 
= fiuj=e2.9 ou lam : : : : 
| Energy Balance 
22 P>5 mm Fiujsk0e y 0.00 10.1 38.9 15.0 57.8 
Recorded 
82 — Si wune|% = 9.9 38.1 13.5 52.0 
Recorded 
60 P<5 mm Socom 14 = 9.7 37.3 12.8 49.3 
22 P>5 mm Recorded = 10.0 38.5 15.5 59.7 


Snowmelt 


The data in category II (P > 5 mm) yield no model fit, which 
suggests that the deviation between recorded and computed snow- 
melt is due to errors in the snowmelt data during heavy rain. 
Austlihylla is a rather flat mountain plateau, but on a larger 
scale the terrain is sloping some 45% (Fig. 2.2). During winter 

a deep snow-drift builds up along the plateau. Water movements 
paraltel to the snow surface may occur, especially during heavy 
rain. These movements, however, have a random character depending 
on physical properties of the snowpack, e.g. ice layers and other 
density barriers, as well as on the rainfall and snowmelt rates 
(Wankiewicz, 1978). The vertical movements of the liquid water 
through the snow will have a corresponding random character. It 
should, however, also be noted that the standard error in the rain 
records is high (Chapt. 3.1). Nevertheless, the conclusion must 
be that during rain-on-snow events a net flow of liquid water has 
passed across the lysimeter edges. Unfortunately, no measurements 
of water equivalent on the lysimeter for discharge verification 


exist. 


By using the 60 days of category I (P < 5 mm), a significant model 
improvement was achieved (Table 4.8). The wind function then found 


was 
f(u) = 2.9 u (wm 2@K7!) 5 (4.2.28) 


where u is the wind speed (ms~!), 


These 60 days largely coincide with a group of days with small 
cloud amounts. During such days the energy balance model at 
Dyrdalsvatn yields higher values than recorded snowmelt, indicating 
that the turbulent heat supply to the snow surface is overestimated. 
It is therefore to be expected that the a-value of the wind: func-~ 
tion, au , during fine weather at Austlihylla is lower than the 
value 3.65 valid during average weather conditions at Dyrdalsvatn. 


Conclusions regarding the wind function. 


The linear wind function (4.2;25) 
flu) = 3.3 u+ 2.3 , 


derived at the open site at Frotveit is not far from optimal at 
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the wooded site as well. Fig. 4.13 illustrates that it is also 

very close to optimal at Dyrdalsvatn (R2=0.79). In addition, con- 
sidering the inaccuracies connected to the deduction of eq. (4.2.27), 
the function (4.2.25) is nearly optimal for the data set at Austli- 
hylla,too. For making the heat flux calculations and drawing the 
figures, the wind function optimal to the respective site is used. 

By exchanging the optimal function with eq. (4.2.25), however, 

very small deviations would occur. The validity of eq. (4.2.25) is 
thus confirmed at several sites, which differ in wind exposure and 
vegetation, and should therefore be suitable in the whole maritime 


region of Western Norway. 


By comparing the wind function (4.2.25) to functions found in 
other investigations, it is essential that the wind speed records 
refer to the same height above the surface. Records at different 


heights can be made comparable by using the formula, 
ni 
7 = A ; (4.2.29) 
where the indices represent two different levels at the same site. 
n varies with such conditions as surface roughness, height, and 
stability. However, the value 0.17 suggested: by Anderson (1976) 


at a melting snow cover is adopted here as an average value. 


In Table 4.9 linear wind functions found in different studies are 
reduced to the reference height, 1.0 m. The adjustable constants 
a and b found in this study are in reasonable agreement with the 
three most reliable investigations (Kuzmin (1961), U.S.A.C.E. 
(1955), Anderson (1976)) listed. 


Table 4.9. Comparison of linear wind function constants found in 


different investigations. 


Value of coefficents for a 
measurement height of 1.0 m 


Reference a(dm 3K" !) b(Wm~¢K~!) 
Kuzmin (1961) based on data from the 3.2 4.0 
USSR plus measurements by Kohler (1950) ie ef 
Snow investigations (U.S.A.C.E., 1955), 4.0 0.0 
Central Snow Laboratory ° ; 
Meiman andGrant (1974), Colorado 1.2. 0.0 
Anderson (1976), using an optimizing 
routine o.4 0.0 
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Scatter diagram of recorded/estimated daily snowmelt 
at Dyrdalsvatn 1979-82 using the energy balance model 
with the optimal wind function f(u) = 3.65u (wm 2k7!) 
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4.16. Daily totals of recorded and modelled snowmelt at 
Dyrdalsvatn during 1979-82. 
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Fig. 4.18. Daily totals of recorded and modelled snowmelt 
at Austlihylla 1980-82. The wind function 


f(u)=3.1 u is used for the energy balance model. 
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The relative importance of the different heat sources in snowmelt. 


The recorded net radiation, Qn ,» and computed sensible (Qi) and 
latent (Q-) heat (using the optimal linear wind function (4.2.24) 
during snowmelt at Dyrdalsvvatn 1979-82 are shown in Fig. 4,18. 
On the average, Qu represented 68%, and Qn 36% of the total 
energy used in melting the snowpack, while Qe represented a 
heat sink of 4%. The data are also separated into two groups 
according to the cloud cover (observed every 3 hour ‘in Bergen) 
(Table 4.10). Cloudy weather often reflects maritime influence, 
while clear or light cloudy weather reflects continental influ- 
ence. Qu ususally provides the highest contribution to the snow- 
melt in both weather types, and is rather independent of the 
cloud cover. Qn is also rather independent of the cloud cover, 
the gain of short-wave radiation during clear weather being » 
counteracted by a longwave radiation loss (Table 4.10). During 
the spring of 1979, when the albedo values were high (74%, Fig.4.1), 
Qu was notably higher than Qn (Fig. 4.19). On the other hand, 
during the last part of the spring of,1980 snowmelt season - 

with high global radiation, Qs , low albedo, a (60%), and low 
wind speed, u - the net radiation, Qn , exceeded the sensible 
heat flux, Qu 


The small size of the latent heat flux, Qe (on the average a 
heat sink of 4% of Quy)» is a consequence of averaging over all 
weather conditions. During 47 days net condensation occurred, 
making a total energy source of 14% of Qy (averaged over the 
122 days). Net evaporation occurred.in 75 days, representing a 
heat sink of 18%. Table 4.10 shows that snowmelt rates on overcast 
days are markedly higher than the rates on less cloudy days. 
Computed energy fluxes suggest that this is mainly due to differ- 
ences in the latent heat flux, Qe . When the weather is cloudy, 
Qp?0 (condensation), while evaporation (Q,<0) dominates in less 
cloudy weather. Evaporation and condensation contribute signifi- 
cantly to the energy budget of the snowpack, while the contribu- 
tions to the mass budget are negligible (Table 4.11). 
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At Frotveit the melting rate was, as a rule, higher at the wooded 
site than at the woodless (Fig. 4.21, Table 4.10). During days 

with light cloud cover (c<0.90), this difference is pronounced, 
while the snow melts somewhat faster at the woodless site during . 
overcast weather (c>0.90). The correlation coefficient, r , between 
recorded daily snowmelt totals at the two sites was 0.82, while 

r = 0.96 using only overcast days (c>0.90), and rv = 0.91 using 
only days when c<0.90. These differences are mainly caused by 
enhanced net radiation, Qn , at the wooded site during clear weath- 
er (Table 4.10)! due to the high contribution of longwave radiation 
from the trees (pp.67). Moreover, the fluxes of Qy and Qe are 
higher at the woodless site than in the wood during humid (over- 


cast) weather. 


The percentage contributions of the fluxes to the total energy 
input at the wooded site are 59% (Qu) >» 32% (Qu); and 9% (Qe) 
while the corresponding values at the open place are 43%, 43%, 
and 14% during the 26 days of simultaneosly recorded snowmelt at 


the two sites. 


Table 4.10 also illustrates that the modelled snowmelt at Frotveit 
and Dyrdalsvatn is overestimated during clear weather, and under- 
estimated during overcast weather. This will not, however, be 
further discussed due to the rather small samples. 
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Table 4.11. Totals of snowmelt, evaporation and condensation (mm 
water equivalent) at Dyrdalsvatn during 1979 and 1980. 


Computed (Net Daily Values added) 


Period Recorded Snowmelt 

Evaporation Condensation 
April |. = 
May 21, 1979 550 mm 10.3 mm 13.4 mm 
April 12, - 
May 1, 1980 291 mm 4.4 mm 1.2 mm 


To conclude, net radiation, Qn , and sensible heat, Qu » on the 
average yielded some 50% of the energy consumed in melting the snow- 
pack in this area, while the contribution from the latent heat flux, 
Qe , was minor. However, the contributions of these three principal 
energy fluxes varied both in space and time. Moreover, it should be 
stressed that the low value of Qe is an average value. During 
shorter periods much heat may be released by condensation or lost by 


evaporation. 
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Fig. 4.19. The variation of the computed energy fluxed, Q 
and Qn , at Dyrdalsvatn during snowmelt 1979-82. 
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Fig. 4.20. Daily mean fluxes of net radiation, Qn » latent heat, Qe ‘ 
and sensible heat, Qu , at the woodless (1), and wooded (II) 
site at Frotveit during the periods of recorded snowmelt 


1981=82. 
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Fig. 4.21. Daily totals of recorded snowmelt at the woodless (1), 
and wooded (I1) site at Frotveit 1981-82. Also indicated 


is the cloud cover, c , observed inBergen. 
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Extreme events 


On May 21, 1979 a snowmelt total of 110 mm was recorded. The melt- 
ing rate was 8 mm hour | in the middle of the day (Fig. 4.22b). 
Fig. 4.22a shows that air temperature, T , and wind speed, u , 
varied throughout the day, but both remained high. The air was moist, 
the albedo low, and the weather cloudy except in the middle of the 
day, yielding a significant contribution from the absorbed global 
radiation, and an effective longwave radiation gain (Table 4.12). 
The rainfall total was markedly lower than the snowmelt total, the 
snowpack in the collection vats was rather small, and the physical 
properties of the snow did not -change on this day. The snowmelt 
measurements should therefore be fairly accurate. When using the 
wind function f,(u) = 3.65 u , the efficiency, R2 , was 0.79 for 
24 hourly observations on that day. Adjusting for the phase lag 
related to the percolation of the meltwater, yields an even better 
fit. Fig. 4.22b suggests a lag of 1 hour, which gives R2=0.88. 


On some days (April 3-4, 1979; Feb. 20-22, 1980; Feb. 3-4, 1982) 

when net radiation loss, rather strong wind, and dry air conditions 
prevailed, no snowmelt was recorded, in spite of diurnal air tempera- 
tures exceeding 0°C. Since the snowpack was approximately ripe and 
jsothermal those days, it can be concluded that no snowmelt occurred. 
The situation may be described as follows: The snow surface gains 
heat through the sensible heat flux, Qu , and gives off heat through 
the net radiation, Qn , and the latent heat flux, Q- (evaporation). 
By suggesting the surface temperature, t , to be OC ils was not 
recorded), the net energy calculated, Qu » is negative, and qs 

must decrease. The energy balance method thus predicts no snowmelt. 
When T, decreases, however, Qy increases, and the net radiation 
loss as well as the evaporation loss decreases, and a balance 
probably occurs (Qy20). The days of February 21-22 may serve as an 
example (Table 4.12). If Te decreases to -1.6°C, a balance is 
established. It should be noted that such a small surface temperature 
decrease only slightly lowers the internal energy of the snow 


because of its poor heat conduction properties. 


These extreme situations illustrate some of the problems in using 
the air temperature as a single snowmelt indicator, and strengthen 
the reliability of the energy balance method. 
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Fig. 4.22a. Recorded wind speed, 
air temperature and 
precipitation at 
Dyrdalsvatn on May 


Zl, U979. 
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Fig. 4.22b. Computed and recorded 
snowmelt at Dyrdalsvatn 


on May 21, 1979. 


qateds 


Snowmelt modelling is of special 
areas for water reservoirs used in 
Meteorological stations are seldom 


Model Sensivity to Data Accuracy. 


interest in outlying catchment 


hydroelectrical power production. 


situated near those areas, and 


the input data to the model therefore have to be extrapolated from 


observations at some distance from 


rest to see how sensitive the model 


thereby introduced. 


Daily mean values of air temperature, T , water vapour, e , 


It is of inte- 


is to errors of the size 


the actual site. 


and 


wind speed, u , are replaced in turns by records from adjacent 


weather stations operated by the Norwegian Meteorological 


Institute. 


T are taken from Bergen-Florida, B (40 m a.s.1.), and Kvamskogen, 


K (408 ma.s.1.), e 
(50 m a.s.1.) (Fig. 
and 


at hs 


Qy 
Florida. 
When comparing the data of x (e, 
Dyrdalsvatn and Bergen-Florida the 


Flesland is chosen as wind station 
The relation between observed wind 


is given by, 


The index 1 


deviations of the data when recorde 
the height corrections are determined without 


method °(2). 
using the data from Dyrdalsvatn. 


second 


from Bergen-Florida, 
The radiation terms 


and u from Flesland, F 
Q, (global radiation), 


(incoming longwave radiation) are taken from Bergen- 


1; Qo, and Q,,) recorded at 
daily mean difference is given by, 


(4.2.30) 


since it is relatively exposed. 
speed at Dyrdalsvatn and Flesland 


(4.2.31) 


indicates the method of determining the systematic 


d at different stations. In a 


and e are thus reduced 


according to the mean height gradients in the area. Fig. 4.23 
illustrates the height gradients in Hordaland, found by using all 
weather stations except the lighthouses, for April 1979. The tempe- 


rature gradient of 


annual mean climatological value in the area (-0.6°C/100 m, 


-0.8°C/100 m is stronger than the accepted 


- O7 - 


Lundquist, 1981). A gradient of -0.8°C/100 m, however, is also 
obtained by using standard normal values (1931-60) for April. The 
reason for the stronger gradient probably is that the more ele- 
vated station sites are snow covered, the lower ones are not. 


The wind speed recorded at Flesland at a height of 10 m is reduced 
to a level of 1.3 m to be comparable with the records from Dyrdals- 
vatn. A surface roughness of 0.15 cm for melting snow (Anderson, 
1976), and a logarithmic wind profile are assumed. It should be 
remarked that this theoretical reduction is rather speculative 
because of the differences in site characteristics, temporal varia- 
tions in surface roughness, and the concept of logarithmic wind 
profile. The rough terrain in the area and the prevailing stable 
stratification above a snow surface make the last concept particu- 


larly questionable. 


Any systematic difference in the cloud cover, andiclear weather: global 
radiation between Bergen and Dyrdalsvatn are neglected. For method 
Ls Qs observed in Bergen is therefore used without corrections 

as model: input at Dyrdalsvatn, and \Qco = 0 . Correspondingly, the 
incoming longwave radiation, Qiy » iS only corrected according to 
the differences in average air temperature by use of Swinbank's 
formula (4.2.12), and 


4Q, = 5.31 + 107 !3(T, - Tp) (wm) (4.2.32) 


Here the air temperature recorded in Bergen, Tp » 18 corrected 
according to the height gradient (Fig. 4.23) to obtain an estimate 
valid for Dyrdalsvatn, Ty 
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4b April 1979 r=- 0.99 
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Fig. 4.23a. Relation between average air temperature of April 1979 
and altitude based on 13 weather stations in Hordaland. 
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Fig. 4.23b. Relation between average vapour pressure of April 1979 
and altitude based on 13 weather stations in Hordaland. 


Table 4.13. Observed (1) and estimated (2) differences in diurnal air 
temperature, AT (°C), humidity, ae (mb), global radiation, 
Q, (wm 2), incoming longwave radiation, AQ y (wm7 2), and 
reduction factor for the wind speed, Au , between Dyrdals- 
vatn and adjacent weather stations at: horizontal distance, 
hd . Standard deviations of the observed differences 


(reduction factor) are given in paranthesis. 


Number of days $; n-= 122 nm = 122 n = 122 n = 122 ni = 53 


Adjacent hd 


weather AT AT Ae Ae Au AU AQ AQ. AQ AQ, we 
stakion k I Z ] 2 ] 2 S So L¥1. Ly¥2 
Bergen- 1 4.0 322 0.83 0.84 0 0 17 17 
Flores (0.9) (0.8) (25) (11) 


Kvamskogen 23 


0.70 0.77 
(0.36) 


Flesland 18 


Method 2 is more general than method 1, since it does not presume any 
measurement of the average difference of the parameters between the 
actual site (Dyrdalsvatn) and the adjacent stations. As can be seen 
from Table 4.13, however, nearly identical corrections are obtained by 
the two methods. The further analysis based on method 1 is therefore 
valid for method 2 as well. (For Kvamskogen, no correction is applied 


in the analysis.) 


Optimal modelled radiation estimates are also used. The global radiation, 
Q. » is modelled according to Table 4.4, where l-c is the cloud para- 
meter. The incoming longwave radiation, Qy » 1S modelled according to 
Table 16, where cf is the cloud parameter, and the albedo, a , is 
estimated according to eq. (4.2.6). A constant value usually proposed, 
a= 0.70 , is also tried. 
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‘Table 4.14a. The efficiency, R2 , of the energy balance model 
when using data from Dyrdalsvatn, D, Bergen-Florida, B, 
Flesland, F, and Kvamskogen, K, for 122 days of recorded 
snowmelt at Dyrdalsvatn 1979-82. Index, M , specifies 
modelled radiation terms. 


Station for recording 


Global Incoming Air Air Wind Average 

Radiation, Albedo poaeeel Temp. humidity Speed time of 
Q. ot Qi y T e u T, @, U R2 
D D B D D D bh 0.79 
D D B D D D 24 h 0.79 
D D B D D F 24 h 0.81 
D D B B D D 24 h O79 
D D B K D D 24 h 0.78 
D D B D B D 24 h 0.83 
D D B D B F 24 h 0.8] 
D D B B B D 24 h 0.84 
D D B B B F 24 h 0.83 
Di Di Diy D D D sh 0.78 
Diy 0.70 Diy D D D ah 0.77 


Table 4.14b. Same as Table 4.14afor 53 days of snowmelt 1980-82 when 


Qry was recorded at Dyrdalen. 
Station for recording i etl 
Us : ay T e u T, @» R2 
: ‘ i Pp D D bh 0.81 
D D D D D D 24 h 0.8] 
D D D D D F 24 h 0.83 
D D D D B D 24 h 0.8] 
: ‘ ‘ B B F. 24h 0.80 
: : ‘ D D D bh 0.73 
: ; : D D D bh 0.77 
B 0.70 B D D D Lh 0.69 
Dn Dy Di D D D bh 0.76 
a ide Di D D D bh 0.70 
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Tables4.14 a and b give computed model efficiencies for different 
input data sets. The first two rows of Table 4.14 a, b differ only 
in the averaging period, 2 h and 24 h, respectively. For both data 
sets R2 is the same. Surprisingly, slightly better model fit is 
found using wu, T, and e-data from stations at a distance of some 
10-20 km from the lysimeter. The difference, however, is not statis- 


tically significant. 


The conclusion must be that the energy balance model is rather 
insensitive to errors in the input data, u, T, and e of the size 
introduced by extrapolating data 10-20 km away when a reasonable 
choise of the wind station is made, and necessary height corrections 
are made. These errors (du = 0.36 u, AT = 0.9°C, se = 0.8 mb, 

Table 4.13) are, however, large compared to the standard errors 

of the data recorded at the actual site (within 0.3 fis", 0.2°C. 

and 0.1 mb, Chapt. 3.1). Since the first-mentioned errors have no 
effect on the model efficiencies, the last-mentioned have none 
either. The accuracy of the measurements of u, 1, and e therefore 


are not critical. 


Table 4,.14shows that the energy balance model is more sensitive 
to extrapolation errors in the radiation input data. The standard 
error of the net radiation, Qn , recorded at Dyrdalsvatn (10 wm? 
Chapt. 3.1) is comparable to the error introduced by substituting 
incoming longwave radiation, Qry , from data recorded in Bergen 
(11-wm2, Table 4.13). Improvements of the radiation data there- 
fore are significant in endeavouring to improve the snowmelt 


estimates. 


Table 4.14aindicates less sensitivity than 4.14 b, probably because 
the data in the former table also include several large snowmelt 
floods when moist and rainy weather invaded the area. In such 


situations, Q is very low and Qy may be rather precisely 


S 
estimated from the air temperature. This implies that during cloudy 
and rainy weather data from the regular network of meteorological 

stations can give fairly accurate snowmelt estimates assuming that 


representative surface wind values can be found. 


~ 102 = 


However, during less cloudy and rainy weather, good estimates of 
net radiation are essential in producing reliable snowmelt estimates. 
Table 4.14 b indicates that the albedo should be estimated by a 


reasonable model if records are missing. 


These results indicate that the rather poor energy balance model 
fit found by Kuusisto (1978) in Finland most certainly was due to 
unrepresentative net radiation records, a point mentioned by the 


author as a possible explanation. 


4.3. The Degree-day Model. 


In the simple degree-day model (3.4.1) the snowmelt is assumed 
proportional to the daily mean air temperature exceeding a threshold 
value. The model is widely used due to its simplicity, and is some- 
times refined by allowing the proportionality factor (degree-day 
factor) to be a function of physical variables. 


For the 122 days of recorded snowmelt at Dyrdalsvatn (1979-82), 
the average value of k in eq. (3.4.1) is 4.0 mm day" ee"", which 
is very close to the value 3.9 mm day lec™!, found by using total 
discharge data from the Dyrdalen precipitation catchment area 
1977-79 (Lundquist, 1980). A rather poor model fit is obtained, 
however (Fig. 4.16, Table 4.15). Table 4.15 also illustrates that 
the degree-day constants varied strongly during the four years of 
snowmelt records. The reason for this instability is the variable 
weather of the area, together with the wind exposure of the site. 
By computing different degree-day constants during cloudy and fine 
weather a considerable improvement is obtained. The degree-day 
factor, k , is very high in cloudy weather, i.e. mild, humid, and 
windy weather effectively melts the snow. The use of average 
k-values therefore markedly underestimates the melting during rain- 


on-snow events during many large snowmelt floods. 
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site } 


Woodless 
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Table 4.16. 
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Degree-day factor, k (mm day 'ec™!), threshold temperature, 
T,(°C), and model fit at Frotveit for the simple degree- 
day method (3.4.1), and the composite model (3.4.2) using 
different values of k and T, when the cloud cover, 
C<0.90 (weather type a), and C>0.90 (b). Also shown is 

the energy balance model fit. 


Model fit 
Model Coefficients Correl.coeff. Model effic. 
k kK. ky Te Tea TH) Be i, r R2 

Degree Day 4.1) 2.6 -—- -—--0.6 *“—- — — — 0.83 0.69 
Degree Day (3.4.2) — 2.3 4.1 -— -0.7 0.5 0.87 0.94 0.90 0.82 
Energy 

By Ane - = = = - = — — 0.85 0.71 
Degree Day (3.4.1) 3.6 — — 0.2 — -= — — 0.93 0.87 
Degree Day (3.4.2) -— 3.8 3.2 — 0.2 0.1 0.96 0.90 0.95 0.90 
Energy 

shee ate - = Se - = — — 0.91 0.83 


Table 4.17. 


“| 
Voc), threshold temperature, 


Degree-day factor, k (mm day 
T,(°C), and model efficiency, R2 , at Austlihylla for the 
simple degree-day method (3.4.1), and the composite model 
(3.4.2) using different values of k and Ty, when the 


rain amount, P<5 mm (category I), and P>5 mm (II). 


Classification Number of days k [? R2 

Simple model (3.4.1) 82 Pat -2.8 0.54 
P<5.0 mm (1) 60 2.3 -2.2 0.73 
P>5.0 mm (IT) 22 Zul fl 0 Gel? 


Composite (3.4.2) 


model 


82 = — 0.58 
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The negative response during February 1980 (Table 4.15) is due to 
the combination of high temperature, very dry air, and net radiation 
loss prevailing in the days of Feb. 21-22 (Table 4.12). 


By correlating daily snowmelt totals recorded at Dyrdalsvatn (1979- 
82) to air temperature, T , observed at Bergen-Florida, a correla- 
tion coefficient of 0.64 is obtained compared to 0.63 when using 

T records from Dyrdalsvatn. The degree-day factor, k , is 

4.0 mm oc lday |, and the threshold temperature, T, , is 3.5°C 
(Dyrdalsvatn: 4.0 mm°C 'day™! and -0.5°C). The systematic deviation 
between JT recorded at Bergen and Dyrdalen (mainly due to the 
height gradient) enters the model through the enhanced threshold 
temperature. A tuned model is rather insensitive to the random 
errors in the temperature field introduced by an intersite distance 
of 10-20 km. The height gradient, however, in general has to be 
estimated. The average deviation between recorded and estimated 
(Chapt. 4.2.3) air temperature difference between Bergen and Dyrdals- 
vatn is 0.8°C. For k = 4.0 mm a this involves an under- 


estimation of 3.2 mm for daily snowmelt totals. This is a signifi- 


day 


cant amount when considering accumulated snowmelt totals, but not 


for flood situations. 


From Table 4.16 (Frotveit) it is seen tnat the model efficiency, 

R2 , at the wooded site is remarkably high, and it is notably lower 
at the woodless site. This can partly be explained by enhanced net 
radiation at the wooded site during fine weather. The contribution 
of the longwave radiation from the trees is closely connected to 

the air temperature. Noteworthy is also the more variable wind speed 
at the woodless site (Fig. 4.2), and the correlation of the wind 
speed difference between the two sites to the cloud cover (humidity) 
(Chapt. 4.1). For these reasons there is a greater model improvement 
at the open site than in the wood when the simple model (3.4.1) is 
replaced with the composite model (3.4.2), and better performance 

of the degree-day method than the energy-balance model at the 

wooded site. However, the data sample is small, making the results 


uncertain. 
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The degree-day factor, k , is highest at the wooded site, mainly 
due to the enhanced net radiation mentioned. Kuusisto (1980) found 
lower values of k in coniferous forests than on open land. This 
illustrates that the constants in general have to be tuned to the 
actual site to be adequate. 


At Austlihylla the performance of the degree-day method is poor 
(Table 4.17). If only the days with precipitation totals of 5 mm 

or less are considered, the days of poorest snowmelt record 

quality are omitted, and the weather variability is reduced. Table 
4.17 clearly shows that the model fit is then improved. (The degree- 
day method performs slightly better than the energy balance model, 
which, however, showed a rather poor model fit. Due to the poor 
quality of the snowmelt records this point will not be pursued.) 


To conclude, the degree-day model may be favourable in regions where 
the climate is stable, and at wind-sheltered sites (wood). The 
optimal degree-day constants, however, show spatial and temporal 
variations, usually necessitating a tuning of the model. 


For wind-exposed sites and areas where the weather during snowmelt 
frequently changes between dry and moist, the simple degree-day 

model will not be suitable. It will in particular fail during snow- 
melt floods when moist and cloudy conditions prevail and in winter 
foehn situations when warm, but very dry and clear air is advected. 
However, different constants during rain and fine weather signifi- 


cantly improve the model estimates. 
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5. Summary and concluding remarks. 


Snowmelt is mainly governed by the weather conditions, in particular 
air temperature, air humidity, wina speed, and net radiation at the 
surface. Reliable knowledge of the relation between snowmelt, 
weather and climate are important in planning hydroelectric power 
production, and for predicting snowmelt floods. Simple models which 
predict melting as a function of recorded or forecasted weather 


therefore are of great practical interest. 


Snowmelt rate estimation is frequently needed in outlying areas, 
usually up in the mountains, where regular meteorological observa- 
tions seldom are carried out. It is therefore important to find 
methods which can use data extrapolated from stations some distance 


away. 


In this investigation snowmelt is studied in a rough, mountainous 
area situated in the maximum precipitation zone some 30 km from 
the sea-coast inWestern Norway and some 10 km from Bergen. The 
climate is maritime, and snowmelt episodes occur throughout the 


snow accumulation season. 


During the years 1979-82 measurements were carried out at three 
locations. The sites at Dyrdalsvatn (1979-82), a nearly level field 
(437 ma.s.1.), and Austlihylla (1980-82) at a sloping ground 

(632 m a.s.1l.) were both situated: above the timberline, while at 
the partly wooded area at Frotveit (1981-82, 285 m a.s.1.) one 
station was located on a woodless site, and a second station in a 


grove of deciduous trees. 


At all the four stations air temperature, air humidity, wind speed, 
and runoff from the snowpack were recorded. In addition, precipita- 
tion, incoming shortwave and all-wave radiation, and albedo were 
recorded at Dyrdalsvatn.At the Frotveit stations incoming shortwave 
radiation and net radiation observations were included, as well as 
albedo and precipitation at the woodless station. Snowmelt (=runoff- 


Z 


rain) was found by lysimetry, measuring runoff from 9 m” plots of 


snow cover. 
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The highest peaks in the snowmelt records during 1979-82 occurred 
during heavy rain. The maximum daily snowmelt was 110 mm, when 

strong wind, high temperature and humidity prevailed, and the sky 
was overcast except for a few hours at midday, and the albedo was 


rather low. 


The collected data are used to test two models for snowmelt estima- 
tion: the energy balance model and the degree-day model. 


’ 


In the simplified energy balance model the energy consumed by snow- 
melt, Qu » is supplied by net radiation, Qn » sensible heat flux, 
Qy » and latent heat flux, Qe 


The net radiation can be directly measured, or it can be determined 


by recording/estimating the single terms of eq. (4.2.1): 


Qy = Qc (l-a) + Ory * eto" : 


where the snow surface temperature, To = 273°K during snowmelt. 


When local topography (shading and surface slope) is taken into 
account, the incoming shortwave radiation may be extrapolated with 
satisfactory accuracy from records from stations some distance away 
(11 km in this study). When such records are not available, eq. 
4.2.44), evolved by using Q< data from Dyrdalsvatn and the fractional 
cloud cover, C , from Bergen, may be applied at sites within the 
maritime region of Western Norway: 


Q, = (-0.16(1-C) + 0.81(1-c)* + 0.071Q,, 


The albedo, ao , during snowmelt is correspondingly modelled by 
eq. (4.2.6): 


@ = -0.13(1-C) - 0.05 In t + 0.87 


Here, t is the number of days which the snow at the surface has been 
exposed to the atmosphere (eq. 4.2.5). The multiple correlation 
coefficient is 0.75. When applied at other climatic regions or at 

at areas where the dust or pollution deposit on the snow is different 
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from that of Dyrdalen (e.g. below the timberline) a tuning of the 


model is required. 


Also the incoming longwave radiation, Qa may be estimated from 


records (R) from stations some distance away when first order cor~ 
rections due to possible systematic differences in the air tempera= 
ture are made, eq. (4.2.32) (Swinbank, 1963): 


-13 776 2 
AQ, = 5.31 + 19713 (T2,-78) (wm 2) 


When reliable records are missing, Qy can be modelled by eq. 
(4.2.115M4): 
2) 


3 


Q, = 1.02 ots + FIC > 92 twine 


which is derived using 150 monthly values of Qry » C and Ti. from 
Bergen. The multiple correlation coefficient is 0.95. The equation 
yields acloud cover term of 71 Wm” 2 For completely overcast weather, 
which is in good accordance with the value 60 Wm 2 given by Palt- 
ridge & Platt (1976, pp. 140). The clear sky term, 1.020T3-92 

(wm 2), yielding 220 wm 2 for T, = 273°K is well in agreement with 
the 229 Wm @ predicted by the equation of Swinbank (1963). -Eq. 
(4.2.11,M4)is probably valid within the whole maritime region of 


Western Norway. 


The net radiation flux should be modelled (recorded) in all weather 
situations except perhaps for snowmelt floods during heavy overcast 
conditions when the term js known for certain to be small. 


The turbulent fluxes of sensible (Qu) and latent (Qe) heat are 
modelled as formulae composed of a wind function factor and a tempe- 


rature/vapour pressure difference between the air and the snow 
surface. Three wind functions were tested against the data set from 
Dyrdalsvatn. The wind functions were optimized by minimizing the 
residual error, o, > between recorded and estimated snowmelt. The 
energy balance model with a theoretically based wind function 
including a stability parameter and a critical Richardson number 


of 0.40 yielded an average surface roughness, Z,of 0.08 cm, which 
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is a reasonable value during snowmelt. The model efficiency, R2 , 
b 


was 0.79. The more simple power law, f(u) = au” , also yielded 
R2 = 0.79 , and the linear wind function, f(u) = aut+b , gave 
R2 = 0.80 . This indicates that the -above-mentioned theoretical 


sophistication does not improve the method, a result whicn is 
probably caused by the deviation from ideal site conditions (surface 


homogeneity). 


The optimal values of a and b in the linear wind function varied 
somewhat in time and space (from one station to another). However, 


during snowmelt the equations, 


Quy = (3.1 U + 2.3)(Ta - To) , (Wm 


1.7(3.1U + 2.3)(e, " e,) >» (Wm 


p=) 
™m 
Ul 


are fairly accurate. The required input dataaredaily values of the 
wind speed, u (ms”!), air temperature, Ta (°C), and vapour pressure, 
e. (mb), measured 1.3 m above the snow surface. During snowmelt the 
surface values are set at To = 0°C , and eg= 6.11 mb. Records 

made within a horizontal distance of some 30 km are adequate when 
first order corrections of T and e (height gradients) are made, 
and u is recorded at a representative wind station. Wind records, 


U made at a level Z, differing from 1.3 m (Z,) should be 


R ? R 

adjusted according to the formula (4.2.29) (Anderson, 1976), 
va. (fay0s17 
ee 


to yield the wind speed values, U at the level Ze 


a b) 
The values of (a,b) in the linear wind function are in reasonable 
agreement with results found by other investigators (Table 4.9). 
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In the degree-day model snowmelt is calculated as a function of -the 


air temperature, T . The simple model,;(3.4.1), 


S. = (mm day!) ; 


m 


koGT * Tal « T 2 Te 
0 , 1* Ty 


where = is the snowmelt rate, is frequently used in snowmelt 
prediction. A composite model also tested is (3.4.2), 


k (7 T seg) » T > Tyq 
during dry (less cloudy) weather 
0 2 J Ty 
*m “Tky(T = Ta)» T > Trp . 
0 » 72 oo during humid (overcast) weather 


In this field study the values of the degree-day constants 
k (mm day !ec7!) and Ty (eco!) at Dyrdalsvatn were 4.0 and 0.5, 
respectively, while the values of (k, ; Tx) were (2.7 , -0.7), 

and (ks , Tao) were (9.4 , 0.5). At Frotveit the values of (k , T,) 
were (2.6 , -0.6) (woodless site), and (3.6 , 0.2) (wooded site). 
These variations illustrate that the degree-day coefficients vary 
considerably with weather conditions and terrain, and that the model 
has to be optimized to the current conditions. The air temperature 
recorded within a horizontal radius of 30 km may give satisfactory 
snowmelt estimates. However, it is important that any height cor- 


rections are precise. 


At Dyrdalsvatn the energy balance model yielded a significantly 
better fit than the simple degree-day method (3.4.1), the model 
efficiency, R2 , being 0.80 and 0.41, respectively. This is due to 
the fact that the energy balance method models the physics of the 
snowmelt, and therefore simulates extreme situations better than 


the degree-day method. 


A significant improvement of the degree-day method at Dyrdalsvatn 
was achieved (R2 = 0.67) by using different degree-day factors 


during overcast and less cloudy weather. 
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At Frotveit the degree-day method yields better fit than the energy 
balance mthod. The variations in the weather conditions (particu- 
larly in the wind speed) were less here, especially at the wooded 
site, and the longwave radiation from the trees is strongly tempera- 
ture-dependent. The correlation between the air temperature and 

the energy fluxes to the surface is therefore high. 


To conclude, the energy balance model is superior when calculating 
snowmelt peaks in maritime regions, and at wind-exposed sites. The 
degree-day method (3.4.1) may give satisfactory estimates in con- 
tinental climatic regions, while it should be replaced by the com- 
posite method (3.4.2) using different degree-day factors during 
humid (overcast) and dry (less cloudy) weather, in maritime climatic 


regions. 


As an average in this study the net radiation, Qn » and the turbu- 
lent flux of sensible heat flux, Qu » each provided about 50% of 

the energy consumed in melting the snowpack, while the contribution 
from the latent heat flux, Qe » was minor. However, the contributions 
of these three principal energy fluxes varied both in space and time. 
On the average, Quy constituted 68%, Qn 36 % and Qe -5% at Dyr- 
dalsvatn, while at Austlihylla Qu and Qn constituted 50% each. 
At Frotveit the percentages were 43, 43, and 14 (woodless site, FI), 
and 32, 59, 9 (wooded site, FII). Except for FI and FII, the periods 
of recorded snowmelt did not coincide. Therefore the contributions 
at Frotveit-Dyrdalsvatn-Austlihylla are not quite comparable. More- 
over, it should be stressed that the low value of Qe is an average 
value. On single days fairly large contributions, mainly positive 
during cloudy weather (condensation) and negative during fine 
weather (evaporation) occurred. The contribution of condensation/ 
evaporation to the mass budget, however, was negligible. 


The distribution of the energy fluxes at Frotveit clearly demon- 
strates the effect of the deciduous wood on the energy balance at 
the snow surface. The canopy porousity to downward radiation flux 
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at the wooded site was estimated to 0.53, corresponding to a sky 
obscuration factor of 0.47. The net radiation is increased by the 
presence of trees, mainly due to longwave radiation, Qy » from 

the canopy. The increase in Qry is greater than the combined effect 
of the decrease in absorbed solar radiation and the slightly lower 
sensible and latent heat fluxes caused by the lower wind speed. 

The difference occurred mainly in fine weather; during cloudy 
weather the net energy input was slightly higher at the open site. 
This is consistent with the observed snowmelt rates, which were 
highest at the wooded site, the difference again occurring in fine 


weather. 


Further research on the topics discussed in this study should 

include .testing/optimizing the wind function constants at other 

sites in different climatic regions. More work should be done regard- 
ing the use of extrapolated/estimated data of air temperature, 
humidity, wind speed, radiation fluxes and albedo. 


When the simplified energy balance model is applied to climatic 

input data, snowmelt flood events may be simulated, and, by a statis- 
tical analysis, the risk for such floods estimated. Attention should, 
however, be drawn to the fact that the science of snow hydrology 
includes many topics not discussed here, e.g. snowpack ripening, 
snowmelt at partly snow-covered ground conditions, water movements 

in snow, ice damming, all of which should be considered when pre- 
dicting areal runoff from (partly) snow-covered fields. 
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